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1.  INTRODUCTION 

This  report  deals  with  some  calibration  problems  of  monitoring  a 
multiple  array  underwater  tracking  range.  The  arrays  in  the  system  are  of  the 
short  baseline  type;  each  contains  four  sonar  transducers  placed  rigidly  at  the 
corners  of  a  cube  in  a  manner  that  describes  a  Cartesian  coordinate  system  in 
three  dimensions.  Figure  1  contains  a  diagram  showing  the  structure  for 
these  arrays  and  indicating  the  numerous  signals  they  may  receive.  It  is 
slightly  deceptive  in  that  real  ray  paths  are  not  straight  lines. 

An  array  receives  a  distinctive  signal  from  a  synchronously  timed  pinger 
attached  to  the  target  vehicle.  The  differentials  of  the  sound  wavefront's 
times  of  arrival  at  the  four  hydrophones  allow  the  computation  of  the 
azimuth  and  elevation  angles  (spherical  coordinate  longitude  and  latitude)  of 
the  normal  to  the  wavefront  at  the  origin  of  the  local  coordinate  system. 
Then,  assuming  direct  path  propagation,  one  can  ray  trace  using  Snell's  law, 
[1],  starting  with  the  aforementioned  elevation  angle  and  utilizing  a  velocity- 
versus-depth  profile  for  the  speed  of  sound  in  the  water.  Finally,  the  time 
differential  between  the  source  pulse  at  the  target  vehicle  and  its  arrival  at  the 
array  is  used  to  stop  the  ray-tracing  algorithm  and  determine  the  location  of 
the  target  relative  to  the  array.  The  local  track  is  the  sequential  set  of  these 
estimated  positions. 

Each  array  in  the  system  operates  over  a  limited  radius.  As  the  target 
sojourns  through  the  range,  it  is  tracked  by  a  number  of  these  arrays.  See 
Figure  2  for  a  plan  view  of  the  Nanoose  range.  (The  zero  level  in  the  vertical 
is  taken  as  mean  sea  level.)  The  overall  path  is  constructed  by  transforming 
each  piece  of  local  track  to  the  coordinates  of  the  range  based  upon  the 
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assumed  location  and  orientation  of  the  various  local  coordinate  systems. 
Discontinuities,  or  mismatches  occur  because  the  track  produced  by  one  array 
does  not  mesh  well  with  that  produced  by  a  neighboring  array  in  the  overlap 
regions.  Such  mismatches  can  be  seen  in  Appendix  B. 

It  appears  that  there  are  several  sources  of  systematic  error  in  the 
operation  of  this  system.  The  question  of  individual  array  location  and 
orientation  has  been  treated  previously  [2].  The  present  report  addresses  the 
timing  synchronization  problem,  from  a  statistical  point  of  view.  That  is,  the 
pulses  received  at  the  arrays  are  timed  to  the  range  clock  with  great  precision. 
The  timing  information  must  be  transmitted  to  the  pinger  prior  to  the  release 
of  the  target  vehicle.  Although  there  is  no  engineering  reason  to  suspect 
noticeable  error  in  this  transfer,  some  data  exhibit  behavior  consonant  with 
such  an  interpretation.  Thus  we  build  a  mathematical  model  to  account  for 
such  a  source  of  systematic  error,  and  develop  statistical  methodology  for 
interpreting  the  data  in  the  light  of  the  model.  Indeed,  if  the  method 
provided  a  way  to  eliminate  mismatches  for  an  entire  run  for  a  single 
vehicle,  there  would  be  great  temptation  to  use  it  as  a  smoothing  filter. 

Generally  there  are  several  sequences  of  time  points,  called  point  count 
sets,  for  which  two  arrays  simultaneously  produce  track.  These  occur  for 
tracking  in  the  overlap  regions  indicated  in  Figure  2.  The  paired  tracking  data 
of  these  point  counts  are  called  "crossover  data."  It  is  the  crossover  data 
generated  as  a  result  of  the  target  vehicle's  entire  trip  that  provide  the 
evidence  suggesting  timing  synchronization  problems  and  the  data  base  for 
evaluating  the  use  of  such  a  model.  Each  version  of  track  in  a  crossover  data 
set  is  assumed  to  have  been  converted  to  the  (common)  range  coordinate 
system.  The  components  of  this  system  are  called  "downrange,  crossrange, 
vertical,"  or  sometimes  "centerline,  crossline,  vertical." 
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Figure  1.  Short  Bose  Line  Arroy  and  Ray  Paths. 
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SCALE  IN  YDS  1x1000) 


The  model  is  developed  in  Section  2.  It  allows  for  both  a  timing  offset 
and  a  drift.  Estimation  methods  and  statistical  properties  are  developed. 
Section  3  contains  applications  to  real  data.  It  is  shown  that  the  offsets  are 
significant  and  the  drifts  are  negligible.  The  issue  of  the  reality  of  the  effects  is 
also  treated.  I.e.,  is  it  the  timing  offset  constant  for  all  crossover  regions  used 
in  a  given  run?  To  answer  this  question,  some  modeling  of  the  components 
of  variance  is  required  and  when  done,  a  statistical  analysis  is  performed.  It 
appears  that  the  true  cause  of  these  effects  is  something  other  than  a  timing 
synchronization  bias.  The  analyses  are  supported  by  some  characteristics  of 
the  noise  process,  and  these  are  presented  in  Section  4.  Conclusions  are 
summarized  in  Section  5.  A  number  of  appendices  are  included  that  contain 
supporting  data  and  information,  including  source  codes  for  computations. 

2.  MODEL  DEVELOPMENT  AND  STATISTICAL  PROPERTIES 

Figure  3  contains  a  mockup  illustrating  conditions  that  support  the 
consideration  of  a  timing  synchronization  correction.  It  shows  a  plan  view 
which  includes  the  radial  lines  from  the  arrays  to  the  estimated  track  points 
in  adjacent  overlap  regions.  Between  these  regions,  much  data  is  supplied 
only  by  array  Aj.  The  analyst  does  not  see  the  radial  lines  on  his  screen,  nor 
does  he  see  the  black  dots.  He  sees  only  the  X's  and  0's  (no  distinction 
between  the  two)  and  no  mismatch  is  apparent.  (Mismatches  would  be 
apparent  however  for  track  pointing  in  a  different  direction.)  But  when  one 
pairs  up  the  radial  lines  by  common  time  points,  then  one  sees  that  the  two 
versions  of  track  lack  coherence  and  can  be  improved  by  stretching  the 
estimated  points  to  the  black  dot  positions.  This  can  be  achieved  by  a  single 
constant  adjustment  to  the  transit  time  values  in  the  ray  tracing  algorithm. 
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Figure  3  :  EffectofTiming  Offset  Errors  (Synthesized) 
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Track  according  to  Arrays  A,  or  A,  as  indicated 


The  figure  riso  helps  one  to  imagine  the  distinction  between  the  effects  of 
an  array  location  error  model  and  a  timing  adjustment  model.  If  the  former 
were  applied  to  the  situation  in  Figure  3,  at  least  two  arrays  would  have  to  be 
moved.  But  any  such  decisions  cannot  be  made  in  isolation.  The  ultimate 
goal  is  the  simultaneous  improvement  of  coherence  in  all  overlap  areas  and, 
of  course,  discontinuities  must  not  be  created  at  other  places. 

A  crossover  data  set  is  a  set  of  matched  pairs  of  three  vectors 


X(t)  = 


X^t) 

X2(t) 

X3(t) 


fY,(t)l 
Y(t)  =  \  V2(t) 

lY?(t)J 


t  =  1,  •••,  T 


(1) 


representing  two  versions  of  the  same  vehicle  track  for  a  common  set  of  time 
values  (point  counts),  T  in  number.  The  array  that  produced  X(t)  is  located  at 
a  and  the  one  that  produced  Y(t'  is  located  at  p.  Thus  these  two  versions  of 
track  can  be  represented  in  their  local  coordinate  systems  as 

£(t)  =  X(t)-a  n(t)  =  Y(t)-p  (2) 


Timing  synchronization  error  corrections  may  be  viewed  as  either  stretching 
or  contracting  the  vectors  £(t)  and  r|(t)  by  the  same  (time)  amount.  The  effect 
of  such  corrections  is  not  constant,  but  depends  upon  (i)  the  speed  of  sound  in 
water  at  the  depths  X3( t)  and  Yj(t);  (ii)  the  elevation  angles  of  the  ray  traces  at 
these  times  and  depths. 

The  magnitudes  of  the  time  adjustments  are  small  and  the  effects  can  be 
described  using  first  order  terms 


*(t,g(t))  =  £(t)  +  g(t)a(t) 
t|(t,g(t))  =  rj  (t)  +  g(t)b(t) 


(3) 


where  g(t)  is  a  scalar  function  of  time  representing  the  total  adjustment  for 
timing  offset  (5)  and  drift  (m) 


g(t)  =  8  +  m(t-to) 


(4) 


for  a  conveniently  chosen  central  time,  to;  and  a(t),  b(t)  are  the  scaled 
directions  of  stretch 


a(t)  =  v(t) 


cos(8(t))cos(<j>(t)) 

<  cos(8(t))sin(<f>(t))  " 
^sin(6(t)) 


(5) 


i.e.,  v(t)  is  the  speed  of  sound  at  depth  ^(t),  0(t)  is  the  elevation  angle  of  the 
ray  trace  from  a  at  ^(t),  and  0(t)  is  the  azimuth  of  £(t)  from  a.  The  vector  b(t) 
is  defined  similarly  for  the  track  q(t)  and  its  origin  p. 

It  is  easily  seen  that,  for  given  5,  m  and  to,  the  corrected  versions  of  track 
in  the  (common)  range  coordinate  system  are 

X(t,g(t))  =  X(t)  +  g(t)a(t) 

Y(t,g(t»  =  Y(t)  +  g(t)b(t) 

The  sta t: stical  estimation  problem  is  to  chose  5  and  m  so  that  the  two 
versions  of  track  agree  as  well  as  possible.  The  least  squares  approach  is 
adopted.  Using  the  notation 

T 

Ave||W(f)f  =-X*V'(/)W(f)  (7) 

1  T  1=1 

we  set  up  the  objective  function 
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(8) 


C  =  Ave|x(i,x(/))-y((,x(f)|2 
=  Ave|ix(( )  -  Y(l)||2  +  Avejglf  )f||«(  f )  -  b(tf 
+2Ave  g(tP(t)-b(t)i[X(t)-Y(t)] 

and  find  the  values  of  5  and  m  that  minimize  Q.  When  the  two  partial 
derivatives  are  set  equal  to  zero  it  is  seen  that  there  is  convenience  in 
choosing  to  so  that 

Ave(t-t0)||fl(/)-b(0|2  =  0  (9) 

That  is 

h  =  Ave  /  Ave 

This  done,  the  normal  equations  may  be  expressed  as 

SAvefrlO-Wtf  =  -Ave(fl(/)-Ml))'[X(l)-y(0] 

(10) 

m  Ave(f  -  (Q  )3||a(()  -  H')f  =  -  Ave(l  -  l0)[n(J)  [X(f)  -  Y(f)] 

and  one  can  readily  solve  explicitly  for  the  minimizing  values  5  and  m. 
Positive  values  of  6  are  associated  with  stretching,  negative  with  contracting. 

If  we  let  Q  be  the  minimizing  value  of  Q,  it  is  useful  to  establish  the 
decomposition 

Q  =  Q  +  (<5-<5)  A\^a{t)-b{t)f  +  Aye{t-t0p{t)-b{tf  (11) 
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Before  justifying  (11)  it  is  convenient  to  shorten  the  writing:  let  W(/,<5(f)) 
=  X(t,g(t))-Y(t ,g(t));  c(t)  =  a(t)-b(t);g(t)=  5  +  m(t-t0).  Then  the  statement 
says 

Ave||w(f,g(f))||~ 

=  Ave|W(f/g(f))|2  +  (5  -  §f  Ave||c(f)f  +  (m  -  m)2  Ave(f  -  t0)\\c{tf  (12) 

and  the  result  is  justified  by  use  of  the  three  orthogonality  relationships 

/ 

Ave|w(f,0)+  <5c(f)+  m(t-  f0)c(f)j  c(t)  =  0 

Ave[w(f,0)  +  Sc(t )  +  m(t  -  t0)c(t)]  (t  -  t0)c(t )  =  0 

Ave(t  -  toy(t)c(t)  =  0  (13) 

which  in  turn  are  established  using  the  normal  equations  (10),  and  (9). 

The  significance  of  the  offset  and  drift  parameters  can  be  judged  if  we 
develop  the  means  and  variances  of  5,  rh.  Letting  E  denote  the  mathematical 
expectation  operator,  we  begin  with  the  assumption 

E[W(t,g(t)]  =  0  (14) 

which  embraces  the  idea  that  the  two  corrected  versions  of  track  produce 
common  track  without  systematic  error. 

Further,  let 

-j-l  r  -j-l 

K{=Xc'(lM<)  ;  X(t-fo)V(fMO  (15) 

_  t  J  L  t 

and  use  (13)  and  (14)  to  show  that  the  estimators  t>  and  m  are  unbiased. 
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E  (S)  =  -K5£c'(i)E[W(f,0)]  =  KjXc'(  l\S  +  m(t  -  f0)]c(!)  =  S 

t  t 

E(m)  =  -K„X(l  -  10>:'(()E[W((,0)]  =  Km£(l  -  f0)c'(l)[<5  +  m(t  -  «„)>(!)  =  m 

t  t 


To  develop  variances,  we  assume  that  the  positive  lag  covariances  are 
zero;  support  for  this  appears  in  the  section  on  noise  characteristics.  Let  M 
represent  the  (zero  lag)  covariance  matrix.  (See  Appendix  C  for  estimates). 


M  =  E[W(/,g(f))W'(f,g(f)] 

(16) 

Then  one  easily  represents 

var(5)  =  K3Xc'(0Mc(0 

t 

(17) 

var(m)  =  K^(t-t0fc'(t)Mc'(t) 

t 

(18) 

cov(£,m)=  KsKm£(f-fo)c'(OMc(f) 

(19) 

t 


If  M  is  proportional  to  the  identity  matrix  then  this  last  term  is  zero;  the  other 
two  terms  simplify  immensely;  and  a  standard  regression  development  can 
be  used.  But  the  study  of  noise  characteristics  does  not  support  this.  On  the 
other  hand  the  vectors  (c(t))  do  not  change  much  with  t  and  this  has  the 
tendency  to  render  (19)  to  be  small.  The  reason  for  this  stability  is  that 
crossover  data  occurs  only  at  the  greater  distances  from  the  sensing  arrays;  the 
azimuth  and  elevation  angles  and  the  layer  sound  speeds  do  not  change 
much. 
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It  appears  that  the  matrix  M  changes  with  the  crossover  data  set. 
Methodology  for  estimating  it  appears  in  Section  4.  Estimates  of  the  M 
matrices  appear  in  Appendix  C. 

3.  DATA  ANALYSIS 

The  data  consists  of  <5,  m  ,  and  supporting  values  for  69  segments  of 
crossover  track  collected  over  four  separate  days  with  two  (temporally  serial, 
not  concurrent)  target  vehicles  (runs)  per  day.  The  information  is 
summarized  in  Table  1.  Missing  variance  estimates  indicate  either  a  data 
shortage  or  outlier  problems.  In  a  few  cases  the  tracks  were  curved,  and  the 
straight  line  model  is  not  adequate.  The  units  are  milliseconds  for  8,  and 
milliseconds  per  point  count  for  m.  One  millisecond  translates  to  about  4.5  to 
5  feet  of  distance.  The  ratios  of  means  to  standard  deviations  are  used  to 
judge  whether  the  effects  are  significantly  different  from  zero.  Virtually  all  of 
the  offsets  and  some  of  the  drifts  are  significant  although  the  latter  are  not 
strongly  so.  The  "r"  column  contains  the  correlations  between  8  and  m  . 
They  are  insignificant.  A  further  search  for  large  scale  drift  was  made  by 
plotting  <5  against  the  crossover  central  time  to  for  each  run.  They  appear  in 
Figure  4.  If  a  smooth  signal  were  discernible  then  we  would  have  a  way  to 
connect  the  5  values  that  appear  in  each  column  of  Table  2.  But  they  are 
chaotic  and  provide  no  incentive  to  continue  any  concern  about  drift.  It  is 
concluded  that  the  offsets  are  significant  and  the  drifts  are  not.  The  latter  will 
be  dropped  from  further  consideration.  Some  graphical  examples  of  the  effect 
our  timing  corrections  have  been  selected  and  appear  in  Appendix  B. 
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Table  1:  Offsets,  Drifts,  Signal  to  Noise  Ratios 


Dale 

Vehicle 

8 

6/Oj 

m 

m/am 

r 

*0 

T 

3/23/89 

MK  30 

1.37 

.0461 

29  62 

.00168 

.001627 

1  03 

-.00 

405  13 

42 

1  74 

.0431 

40.28 

-  00831 

.002217 

3.75 

.00 

422.72 

33 

1 ,69 

.0925 

18  29 

-.01559 
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Table  2:  Estimates  of  the  Timing  Offset 
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Figure  4.  Offset  vs.  Central  Point  Count  for  the  eight  runs 


A 

Since  the  8  values  are  significant,  it  is  important  to  discover  what  they 
really  represent.  They  have  been  labeled  timing  synchronization  offset  errors 
and,  if  that  labeling  is  a  valid  physical  description,  then  each  run  (vehicle)  on 
each  day  should  have  its  own  unique  value  and  this  common  value  could  be 
used  to  correct  all  track,  not  just  crossover  track.  What  follows  is  an  analysis 
of  this  point. 

Table  2  contains  a  graphical  positioning  of  the  sixty-nine  8  values 
according  to  the  eight  runs.  There  were  eleven  array  pairs  involved  in  this 
process  and  they  are  marked  with  the  letters  a,  b,  ...,  k.  (The  correspondence 
with  the  overlap  regions  marked  in  Figure  2  is  given  in  Table  3.)  Study  of 
Table  2  shows  considerable  "within  run"  variability,  and  one  is  inclined  to 
doubt  the  reality  of  the  constant  offset  interpretation. 

TABLE  3.  IDENTIFICATION  OF  ARRAY  PAIRS  TO  RECONCILE  TABLE  2 


WITH  FIGURE  2 


a 

(1,11) 

e 

(3,4) 

i  (5,56) 

b 

(1,2) 

f 

(43) 

j  (5,14) 

c 

(2,11) 

g 

(5,6) 

k  (15,16) 

d 

(2,3) 

h 

(435) 

In  passing  we  also  note  that  the  8  values  for  each  array  pair  appear  to 
have  some  temporal  coherence,  and  hence  the  assignable  causes  may  be 
related  to  this,  but  that  is  an  issue  for  another  time. 


We  proceed  to  model  the  components  of  random  error  variance  and 
develop  a  test  statistic  for  examining  whether  or  not  all  8  for  the  same  run 
can  be  viewed  as  constant.  The  {<5, }  are  modeled  as  being  affected  by  the  day 
(water  depth  velocity  profiles  change  from  day  to  day),  the  run  (the  second 
run  is  a  different  vehicle  operating  later  in  the  day),  and  the  array  pairs 
generating  the  crossover  data.  The  array  pairs  are  assumed  to  produce  fixed 
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effects,  with  values  <Xj  for  j=l,  ...,  11;  and  the  day  and  run  effects  are  assumed 

2  2 

to  be  random  with  zero  means  and  variances  aj  and  o2  respectively.  There  is 

2 

also  an  experimental  (residual)  error  having  zero  mean  and  variance  o  . 

The  result  is  the  mixed  model 

8  =  Xa  +  ZP  +  e  (20) 

where  5  is  the  69  component  column  vector  of  6  estimates;  X  is  a  69  x  11 
matrix  of  zeros  and  ones  (incidence  matrix)  relating  the  <5  component  to  the 
array  pair;  a  is  the  11  component  vector  of  array  pair  fixed  effects;  2  is  a  69 
rowed  partitioned  incidence  matrix 

Z  =  (Z1/Z2)  (21) 

with  Z\  having  4  columns  relating  the  8  component  to  its  day  and  Z2  having 

.  >  > 

8  columns  relating  the  8  component  to  its  run.  The  vector  P'  =  (Pr  p2) 

represents  the  random  variables 

'  2 

Pi  =  (Pn,  Ph)  independent  Ntt),©]) 

'  2 

P2  =  (P21/  •••/  P28)  independent  N(0,g2) 

and  the  residuals  are 

e'  =  (ei, e69)  independent  N(0,o2) 

Further  the  vectors  Pi,  p2;  e  are  assumed  to  be  mutually  independent. 

We  will  be  applying  the  model  (20)  separately  to  each  of  the  eight  runs  for 

purposes  of  estimating  the  fixed  effects  (a)  for  each  of  those  runs.  Such 

2  2 

estimates  will  require  values  for  the  variance  components  (o2,  0],o2).  We 
prefer  to  estimate  them  but  once,  by  pooling  all  of  the  data.  This  can  be 
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accomplished  separately  from  the  estimation  of  fixed  effects  by  using  the 
restricted  maximum  likelihood  method:  the  unequal  block  size  version  has 
been  treated  by  Patterson  and  Thompson,  Ref.  [4]  and  is  used  here.  (See 
Appendix  A  for  a  description  of  the  algorithm.)  The  results  are 

o?  =  0.3571  02  =  0.7207  o2  =  0.2430  (22) 

and  the  units  are  seconds  squared.  It  is  noted  that  the  run- to-run  variance  is 
double  that  of  the  day-to-day  variance,  which  in  turn  is  larger  than  the  error 
variance  by  a  ratio  of  about  seven  to  five.  But  the  degrees  of  freedom  for  the 
first  two  variances  are  small  and  the  estimates  may  not  be  very  reliable. 

Using  the  model  (20)  one  sees  that  the  covariance  matrix  of  the 
observables  is 

H  =  c2In  +  ZTZ'  (23) 

where  In  is  the  identity  matrix  of  order  n,  T  is  diagonal,  the  first  four  values 
2  2 

being  O]  and  the  next  eight  02;  and  Z  is  partitioned  (Zj,  Zj)  as  before.  For  the 

full  data  set  n  =  69,  but  recall  that  our  goal  is  to  check  whether  the  fixed  effects 
(offset  estimates)  can  be  viewed  as  constant  within  each  run.  The  plan  is  to 
estimate  the  fixed  effects  for  each  run  and  test  whether  thev  can  be  viewed  as 
constant  (over  the  array  pairs). 

To  do  this  we  proceed  as  follows.  The  6  values  for  the  kth  run  can  be 
identified  by  the  ones  in  the  kth  column  of  Z2,  k  =  1,  ...,  8.  The  array  pairs 
involved  (not  necessarily  all  eleven)  are  identified  by  ones  in  the  X  matrix 
restricted  to  the  k,h  run.  The  Z  matrix  is  restricted  in  a  like  fashion  and  the 
covariance  matrix  (of  order  n,  the  number  of  crossover  data  sets  in  the  k,h 
run)  has  the  same  structure  as  (23)  with  the  reinterpretation  of  inputs.  Of 
course  the  estimates  (22)  must  also  be  input. 
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Let  K  be  the  number  of  array  pairs  involved  in  the  kth  run.  Use  the 
Aitken  estimator  [3] 

dk=(x'kHi'xk)~'x‘kHi'Sk  (24) 

and  its  covariance  matrix 

varfd^JxiHt’XtJ'Uv-'  (25) 

and  the  subscript  k  modifies  the  previous  definitions  of  symbols  so  that  only 
the  kth  run  is  involved.  If  the  null  hypothesis  is  true  (i.e.,  the  ctk  all  fall  on 
the  main  diagonal  of  the  K  dimensional  space),  then  the  maximum 
likelihood  estimate  of  the  common  value  is 

a  (ze>) 

>'= ij=i  i  i 

It  also  follows  from  the  normality  assumptions  that  the  quadratic 

(ak -a)  V(ak -a)  ~  Chi  Square  (K- 1)  (27) 

and  this  statistic  is  the  weighted  distance  of  the  components  of  &  from  the 
main  diagonal  (i.e.,  the  constant  fixed  effect  that  is  the  same  for  all  array 
pairs).  Thus  the  null  hypothesis  should  be  rejected  when  this  distance  is  too 
large. 

The  numerical  results  for  the  eight  runs  are  contained  in  Table  4.  The 
column  marked  p*  contains  the  p  values  (empirical  significance  level) 
corresponding  to  the  test  statistics  listed  under  the  "distance"  column.  Seven 
of  the  eight  values  indicate  rather  rare  events  and  the  eighth,  p*  =  0.14,  is 
associated  with  one  degree  of  freedom.  Low  degrees  of  freedom  tend  to 
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dampen  the  chances  of  finding  significant  results.  The  final  column  is  the 
ratio  of  the  distance  (27)  to  its  standard  deviation  under  the  null  hypothesis. 

We  assert  that  the  array  pair  effects  are  not  constant  and  that  there  are 
other  sources  of  systematic  error  that  dominate  this  process. 


TABLE  4.  TESTING  THE  EIGHT  RUNS  FOR  CONSTANT  FIXED  EFFECTS 


a 

K-l 

Distance 

P* 

; 

distance/SD 

0.61 

6 

30.17 

3.6xl0-5 

8.71 

3.03 

1 

2.12 

0.14 

1.51 

2.35 

8 

17.71 

0.06 

4.43 

3.09 

9 

51.67 

5.2x1 0-8 

12.18 

3.05 

4 

41.48 

2.1x10-8 

14.67 

2.13 

6 

14.16 

0.03 

4.09 

1.79 

7 

18.13 

0.01 

4.85 

0.46 

8 

30.44 

1.8x10-4 

7.61 

4.  CHARACTERISTICS  OF  NOISE 

In  order  to  study  the  stochastic  nature  of  the  tracking  errors  we  fit  straight 
lines  to  the  track.  (This  was  done  even  for  track  that  did  not  appear  linear  via 
a  visual  scan  of  graphical  output.  The  exceptional  cases  are  marked  in 
Appendix  C.)  It  is  assumed  that  the  target  vehicles  speed  is  constant.  The 
result  is  a  set  of  deterministically  spaced  points  that  fall  on  a  straight  line  in 
three  space.  From  these  we  can  produce  residuals  and,  assuming  local 
stationarity,  study  the  noise  structure. 
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The  basic  (unadjusted  for  constant  speed)  straight  line  is  found  by 
principal  components.  Let  X(t)  be  used  for  X(t,0),  eq.  (6)  and 

X=7lX(f)  (28) 

1  l 

We  seek  a  projection  (direction)  p  =  (pi,  Pb  P3)' 

Z«-5>/Xi(0  =  P'X(t)  (29) 

1 

that  will  maximize  the  variance  of  (Z(t)} 

CTz  =  7l  [m-z]2  =  p'c,p  (30) 

1  t 

where 

(X(/)-X)(X(f)-X)  (31) 

‘  t 

is  the  covariance  matrix  of  the  track  data,  X(t). 

Since  the  vector  p  is  a  set  of  direction  numbers,  we  adopt  the  usual 
constraint,  ^pf  =  1-  (This  will  keep  o \  finite.)  It  is  easily  shown  that  the 
three  solutions  to  the  eigen  problem 

Cxp  =  A  p 

provide  us  with  an  orthonormal  basis 

P  =  {P,,P2,P3}  (32) 

2 

where  P,  is  the  eigen  vector  corresponding  to  X,  and  X]  >  >  X3.  Since  07  = 

2 

p'C*p  =  Ap'p  =  X,  and  our  goal  is  to  maximize  o^,  we  choose  the  first  eigen 
vector  for  use. 
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The  set  of  values 


Z,(f)=P,'(X(l)-X)  (33) 

represents  the  succession  of  values  after  the  vectors  [X(f)-Xj  are  projected 

onto  the  line  of  the  first  principal  component.  These  projections  are 
orthogonal  projections  and  may  not  well  represent  the  relative  positions  of  a 
target  vehicle  heading  in  the  direction  Pi  at  constant  speed.  Accordingly  it  is 
appropriate  to  make  an  adjustment  in  the  (Zi(t)}  to  account  for  constant 
speed.  The  collection  of  point  counts, 

f]  <  tj  <• .  •<  tj,  (34) 

all  differ  by  integral  multiples  of  some  base  value,  A,  representing  distance 
traversed  per  point  count.  Let  us  perform  a  simple  linear  regression  of  the 
Zi(t)  on  the  times  (34).  The  fitted  values  are 

Z{t)  =  a  +  b(t-i)  (35) 

where  fl  =  —  ^(Z](t)  and  h  =  ]T((f-F) Z3(f) /^((f - f )“.  Now  the  successive 

values  of  Z(t)  dif ter  by  integral  multiples  of  a  converted  base  value,  bA,  and 
represent  the  constant  speed  progression  of  the  vehicle  in  the  direction  P3. 

It  remains  to  represent  the  distances  Z(f)  in  the  original  coordinate 
system.  Let 

R(t)  =  (z(/),0,o)'  (36) 

be  the  estimated  position  of  the  vehicle  in  the  basis  P,  and  then 

Xsj(0  =  X  +  PR(t)  (37) 
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will  be  the  straight  line  in  the  original  coordinates.  Similarly,  when  (Y(t)}  is 
put  into  this  algorithm,  we  produce  Vsj(t). 

Now  we  are  positioned  to  estimate  the  covariance  matrices  of  the  noise 
processes 

Dx  =  fX[X(0-X5,(!)][X(t)-Xs,(!)]- 
1  t 

°y = ^  Iim- (as) 

1  t 

DxJ/=7l[X(')-Xs,(t)][Y(l)-7'sl(0]' 

'  t 

Computational  work  shows  considerable  variability  in  these  matrices  as  one 
changes  the  day,  the  run,  and  the  array  pair.  Also  the  cross  correlations  are 
mostly  different  from  zero.  The  matrix  M,  eq.  (16),  is  estimated  by 

M  =  Dx  +  Dy  -  D^y  -  Dxy  (39) 

and  used  in  (17),  (18),  and  (19).  The  quantities  (38)  and  (39)  appear  in 
Appendix  C,  and  illustrate  their  variable  nature. 

Another  immediate  use  of  the  noise  processes  is  to  look  at  the 
autocorrelations.  In  addition  to  the  three  individual  components  we  define  a 
"noise  displacement"  process 

l 

•*(0  =  ji(x,(<)-Xslj(f))2j2  (40) 

(Of  course  the  same  is  done  for  (Y(t)J.) 

The  autocovariances 
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(41) 


1  t=l 

for  (d(t)}  and  for  the  individual  components  {Xj(t)-XS].(t)}  are  computed  and 

normalized.  Plots  of  R(h)/R(0)  vs  h  appear  in  Appendix  D.  It  is  typical  that 
they  become  and  remain  in  a  general  level  of  "static"  for  h  >  1. 

5.  SUMMARY 

The  paper  contains  methodology  for  estimating  the  presence  of  timing 
synchronization  error  and  judging  its  significance  in  the  framework  of  our 
short  baseline  array  underwater  test  range.  The  methodology  was  applied  to 
real  data.  Some  plots  illustrating  the  effects  appear  in  Appendix  B.  One  must 
view  them  with  care  because  the  coordinate  scales  are  so  different.  It  is  rather 
typical  that  some  unresolved  systematic  error  is  exposed  in  the  side  view. 

Analysis  of  the  results  does  not  support  the  timing  synchronization  error 
model  as  accounting  for  the  discrepancies.  Other  sources  of  systematic  er  ror 
must  be  unmasked  first.  Possibilities  include  the  orientation  of  the  arrays  [2], 
biases  in  the  raytracing  inputs  [3],  and  perhaps  a  temporal  or  spatial  gradient 
in  the  depth  velocity  profile. 

Appendices  C  and  D  contain  information  about  the  second  order 
properties  of  the  three  dimensional  noise  process  and  about  time  lag 
correlations. 
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APPENDIX  A.  PATTERSON  AND  THOMPSON  ALGORITHM 


The  purpose  of  this  appendix  is  to  document  the  computational  method 
for  estimating  variance  components  that  is  but  implicitly  described  in  the 
Patterson  and  Thompson  paper  [4].  We  use  the  notation  of  that  paper. 
Specifically 

y  =  Xa  +  e  (A.l) 

where  the  observeable  y  is  an  n  component  column  vector,  X  is  an  n  by  t 
matrix  of  rank  t  determined  by  the  allocation  of  treatments  to  units,  a  is  a  t 
vector  of  fixed  effects,  and  e  is  a  mean  zero  n  vector  of  normal  random 
variables  with  covariance  matrix 

V  =  o2H  H  =  ZrZ'  + 1.  (A. 2) 

The  matrix  I  is  the  identity  of  order  n,  Z  is  the  n  by  b  design  matrix  for  c  block 
factors,  and  I'  is  a  diagonal  matrix  containing  the  variance  components 
relative  to  the  basic  error  variance  a2.  I.e., 

C 

H  =  I  +  XZpZp'Yp  (A. 3) 

p=i 

and  each  block  design  matrix,  Zp,  is  n  by  bp;  each  plot  having  exactly  one  level 
in  each  Zp,  p  =  1, ...,  c,  the  variance  of  the  p,h  block  is 

2 

Op  =  Ypo2  for  p  =  1,  2 c  (A. 4) 

and  Z  in  (A. 2)  has  the  partitioned  form 

Z  =  (Z,,Z2/...,ZC).  (A. 5) 
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The  diagonal  elements  of  T  are  b]  consecutive  -ft's,  b2  consecutive  y2's,  and 

C 

bc  consecutive  yc's,  with  Xbp  =  b. 

p=i 

Out  goal  is  the  estimation  of  yi,  y2, yc,  a2  by  the  restricted  maximum 
likelihood  method  for  the  unbalanced  case.  The  steps  follow.  Let 

S  =  I  -  X(X'X)-1X'  (A. 6) 


W  =  Z'SZ+T  1 


(A. 7) 


The  algorithm  is  an  iterative  one,  and  initial  values  for  the  yi,  ...,  yc  are 
needed  to  develop  W  in  (A. 7).  The  b  by  b  matrix  W  and  its  inverse  can  be 
partitioned  according  to  bi,  b2, bc-  So  doing  allows  both 

p  =  W_1Z'Sy  (A. 8) 

u  =  r’-r’w’r1  (a.9) 

to  be  viewed  as  partitioned.  I.e., 


P  -  IP]/  P2’  •••'■Pc 


(A. 10) 


u  =  (Ujj)  i,  j  =  1, ...,  c 


and  Uij  is  a  b,  by  b,  matrix. 

Now  we  are  positioned  to  define 


R  =  y'Sy-y'SZp 

Pp  =  Pp  Pp/Yp 


(A. 11) 


and  Ep  =  trace  (Upp)  for  p  =  1,  ...,  c.  Then  define  the  modified  information 
matrix  {fj,},  of  order  c+1, 


fij  =  trace  (UijUji)  for  i,  j  =  1, c 


fp,c+i  —  fc+i'p  -  Ep  for  p  —  1, c 

fc+l,c+l  =  tt—t 

and  let  {f‘i}  be  its  inverse. 

Finally,  the  estimator  update  equations  are,  using  k  =  c+1 

d2  =  f^f^Bj+f^R 
;=i 


(A. 12) 


7,  =  f!°’  +  -tj-]  XfB,  +  fkR  for  i  =  . . 

^  l)-' 


(A.13) 


and  the  new  {/,)  can  be  placed  in  (A.7)  to  start  the  next  iteration.  No  initial 
value  for  o2  is  required,  but  it  must  be  updated  at  each  iteration.  The 
algorithm  should  be  stopped  when  (A.13)  is  stable.  The  variance  components 
are  computed  from  (A.4) 

The  application  in  the  present  report  has  c=2,  with  bi  =  4  and  b2  =  8.  It  is 
helpful  to  record  the  partitioning  and  inversion  of  W,  (A.7). 


w=  zjSZ,  Z',SZ2  +r/,/7,  o 


z2sz]  z2sz2 1  L° 


J2ly2 


(A.14) 


Wn  vvI2 

W21  w22 


symmetric. 


rwn  vv12 

w1  - 

w21  w22 


(A.15) 


can  be  computed  from 
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w22  =  [w22  -  W^wrfw^]"1 

w12  =  -Wf^W^W22  (A.  16) 

w11  =  [;1-w12w21]wf11 
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APPENDIX  B.  EFFECTS  OF  CORRECTIONS 


Some  selected  results  of  applying  the  method  are  given  in  Figures  B1 
through  B4;  two  solutions  for  each  of  the  four  days,  top  and  side  views  for 
each.  Original  track  from  the  first  member  of  the  array  pair  is  marked  with  a 
cross,  and  from  the  second  member  with  a  small  circle.  The  (timing) 
corrected  tracks  are  marked  with  connected  solid  lines  for  the  former  and 
dashed  lines  for  the  latter. 

The  first  set  in  Figure  B1  shows  no  noticeable  corrections,  8  is  small.  The 
second  set  in  that  figure  has  8  =  3.54  and  is  a  real  data  version  of  the  situation 
depicted  by  Figure  3.  The  first  set  of  Figure  B2  appears  as  if  one  track  is 
corrected  more  than  the  other.  But  this  can  be  explained  because  the 
"stretching  angles"  are  different.  The  second  set  of  this  figure  shows  a  rather 
common  condition  in  that  the  corrected  track  is  at  a  shallower  depth  than  the 
original. 

The  first  set  in  Figure  B3  has  something  of  a  "showcase"  flavor;  the  top 
view  has  corrected  track  with  desirable  coherence.  (But  a  systematic 
separation  in  the  vertical  remains.)  The  second  set  also  shows  tantalizing 
improvement.  The  first  set  in  Figure  B4  has  8  =  -0.8,  and  provides  an 
example  of  modest  "contracting"  of  points.  The  second  set  is  an  extreme  case 
of  curved  track.  The  top  view  has  an  excellent  correction  displayed,  but  the 
side  shows  that  the  vertical  is  still  separated. 
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APPENDIX  C.  COVARIANCE  MATRICES 


The  covariances  of  the  residual  processes,  eq.  (38),  are  recorded  here  along 
with  M,  eq  (39),  which  estimates  M,  eq.  (16).  The  cross  covariances  Dxy  are 
converted  to  correlations,  Rxy,  by 

l 

Kay(U)  =  DxV{i>j)/[Dx(i’')Dy{j>i)]2-<  W*  1,2,3. 

The  order  of  presentation  is  that  of  Table  1.  The  number  of  points  used,  T,  is 
occasionally  smaller  than  its  Table  1  counterpart.  This  is  due  to  outlier 
rejection  based  upon  a  visual  view  of  the  track.  In  three  instances  these 
values  are  zero  and  then,  the  covariances  are  not  computed. 

In  a  number  of  instances  the  covariances  are  marked  with  a  "c"  or  "cc"  in 
the  identifying  column.  This  means  the  track  was  curved  (c)  or  excessively 
cu-ved  (cc)  so  that  the  straight  line  filter  could  not  reasonably  be  assumed  tc 
provide  a  valid  representation  of  the  path.  It  is  a  curiosity  that,  in  many  of 
these  instances,  the  cross  correlations  are  strong  and  their  use  produces,  via 
eq.  (16),  compensatory  values.  That  is,  the  M  matrices  do  not  appear  to  be 
especially  large.  Some  large  cross  correlations  also  appear  in  cases  that  pass 
the  visual  test  for  straight  line  track. 
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Table 

C.  1. 

Covariances  and  Cross 

Correlations 

of  Residuals 

Date 

T 

Dx 

Dy 

Rxy 

M 

3/23/89 

42 

r 

MK-30 

12.91 

1.46 

-5.27 

5.47 

4.01 

-4.69 

.39 

.85 

-.12 

11.84 

2.65 

-7.93 

1.46 

1.24 

-5.06 

4.01 

6.43 

-.46 

.07 

.12 

-.08 

2.65 

6.97 

-3.88 

-5.27 

-5.06 

23.17 

-4.69 

-.46 

22.77 

-.04 

-.13 

.21 

-7.93 

-3.88 

36.51 

33 

4.26 

-4.82 

6.59 

4.06 

4.42 

-7.29 

.47 

.31 

-.49 

4.43 

-.50 

.42 

-4.82 

7.50 

-4.08 

4.42 

7.06 

-3.50 

-.29 

-.25 

.23 

-.50 

18.26 

-12.38 

6.59 

-4.08 

16.36 

-7.29 

-3.50 

23.17 

.37 

.18 

-.48 

-.42 

-12.38 

58.27 

19 

r 

r 

14.29 

-24.55 

13.39 

I  2.23 

4  SO 

-1.62 

.26 

.48 

.12 

13.55 

-22.53 

9.19 

-24.55 

49.68 

-11.93 

4.60 

14.02 

4.31 

-.17 

-.33 

-.16 

22.53 

81.28 

-8.68 

13.39 

•11.93 

29.00 

-1.62 

4.31 

14.98 

.11 

.16 

-.10  , 

9.19 

-8.68 

47.95 

45 

■ 

r 

■ 

11.73 

.49 

-3.60 

4.17 

-5.00 

4.84 

.32 

-.62 

.53 

11.46 

-3.24 

-2.48 

.49 

.71 

-2.89 

-5.00 

7.78 

-2.83 

.02 

.01 

.09 

-3.24 

8.47 

-6.21 

-3.60 

-2.89 

13.00 

4.84 

-2.83 

12.36 

L 

-.01 

-.14 

-.20  j 

-2.48 

-6.21 

30.39 

41 

2.37 

3.36 

-1.75 

16.93 

1.77 

2.34 

■ 

.06 

.01 

.07 

18.56 

6.17 

.09 

3.36 

11.96 

10.21 

1.77 

1.37 

-4.09 

-.61 

-.10 

.'•2 

6.17 

14.13 

3.65 

-1.75 

10.21 

26.59 

2.34 

-4.09 

19.65 

-.15 

-.04 

.11 

.09 

3.65 

41.39 

22 

, 

. 

. 

r 

- 

4.10 

6.89 

3.70 

5.40 

-10.92 

-1.90 

1 

j 

-.23 

.43 

-.26 

i  11.63 

-5.77 

4.49 

6.89 

14.24 

',.59 

-10.92 

28.52 

-7.00 

! 

-.18 

.33 

-.18 

!  -5.77 

29.41 

-5.05 

3.70 

1.59 

1  .71 

-1.90 

-7.00 

19.49 

1 

-.07 

.17 

-.18 

1  4.49 

1 

-5.05 

36.67 

21 

8.60 

r 

-i 

5.17 

•11.92 

!  5.38 

2.38 

-17.76 

.08 

-.00 

-.07  | 

9.70 

-8.10 

-7.91 

C 

-11.92 

30.49 

-18.51 

i  2.38 

14.14 

-1.35 

-.16 

-.02 

.33 

-8.10 

45.39 

-25.45 

8.60 

-18.51 

16.16 

1-17,76 

J 

-1.35 

68.60 

-.03 

.01 

.09 

-7.91 

-25.45 

78.93 

3/23/89 

37 

1 

4.49 

-13.34 

MK-30 

2.56 

-2.78 

5.52 

5.26 

.25 

-.12 

-.48 

1  6.22 

4.57 

-4.08 

-2.78 

6.39 

.62  ! 

5.26 

12.43 

-1.50  | 

-.25 

-.11 

.27 

l  4.57 

20.70 

- 1 .07 

! 

4.49 

.62 

18.37  1 

-13.34 

-1.50 

49.73 

! 

.00 

-.22 

-.22 

,  -4.08 

-1.07 

81.51 

J 

L 

l 

47 

- 

n 

r 

1 

.11 

r 

.92 

-1.52 

1.18 

i  10.33 

17.40 

-14.20 

-  .14 

-.10 

!  12.11 

15.42 

-13.82 

-1.52 

4.53 

.98 

j  17.40 

32.04 

-2C.33 

1 

.21 

.19 

-  .11 

:  15.42 

3'. 92 

-18.64 

1.18 

.98 

6.06 

j - 14 .20 

-20.33 

25.14 

1 

[ 

-.02 

.08 

•°9  j 

|  13.82 

-18.64 

29.0- 

5/10/89 

43 

; 

MK-27 

i 

3.77 

.50 

-.38 

!  2.30 

.96 

-.13 

j 

.31 

-.35 

.02 

4.25 

1 .00 

-  ,ii 

.SC 

.89 

-1.57 

j  .96 

1.82 

.16 

i 

.36 

.03 

.09  ; 

1.00 

2.64 

-2. -3 

j 

-.38 

-1.57 

4.86 

1  ",5 

.16 

1.82 

i 

-.18 

.34 

-.29  j 

-.11 

-2.13 

8.39  ; 

J 

j 

L 

38 

1.14 

-1.25 

1.38 

l 

1.09 

.87 

-1.08 

.53 

.29 

-.36  j 

1.04 

-1.07 

.30  ] 

-1.25 

3.21 

-.86 

.87 

1.46 

-.17 

.12 

-.08 

.n  j 

-1.07 

5.03 

-2.13 

1.38 

-.86 

3.11 

( 

-1.08 

-.17 

2.40 

.40 

.31 

-22  J 

.30 

-2.13 

6.74  | 

J 

J 

i 

J 

25 

i 

2.05 

1 

*1 

„ 

-3.96 

.96 

.48 

.68 

-.45  1 

| 

-.26 

-.31 

.14  i 

;  3.06 

-2.90 

1.04  i 

i 

-3.96 

9.36 

1.44  ' 

'  .68 

3.02 

1.87 ; 

j 

.11 

.26 

.05  ! 

'  -2.90 

9.66 

3.40 

{ 

.96 

1.44 

7.37  j 

;  -.45 

1.87 

4.08  1 

-.27 

-.05 

,32j 

-i 

L 

1.04 

3.40 

7.98 

34  - 

2.93 

r 

i 

1 

-.16 

.32 

1.11 

-1.48 

1.31 

j 

-.07 

.36 

.11  • 

4.31 

-1.77 

.59 

i 

-  .16 

.45 

-1.15  | 

1  -1.48 

3.33 

-i.o9  ; 

i 

-  .04 

.43 

-.08 

-1.77 

2.75 

-  .96 

i 

.32 

-1.15 

4.75  1 

1  1.31 

-1.09 

3.13 

.26 

-.79 

.38 

.59 

-  .96 

4.99 

-»  L  J 
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Date 


CC 


5/10/89 

MK-27 


C 


C 


T 

Dx 

Dy 

Rxy 

56 

1 

2.63 

-3.61 

1.94 

1.93 

1.84 

-1.72 

.32 

-.33 

.05 

-3.61 

6.23 

-1.96 

1.84 

5.70 

-2.40 

-.15 

.18 

-.11 

1.94 

-1.96 

4.08 

-1.72 

-2.40 

3.09 

-.01 

-.22 

.01 

59 

_ 

2.29 

-2.16 

3.26 

1.35 

.28 

-.81 

.33 

-.40 

.31 

-2.16 

4.59 

-1.95 

.28 

2.78 

-1.14 

.03 

.21 
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3.26 

-1.95 

7.30 

-.81 

-1.14 

4.27 

-.01 

-.28 

.40 

50 

. 

132.12 

525.44 

-7.44 

109.42 

465.46 

9.12 

.96 

.79 

.59 

525.44 

3906.70 

-7.16 

465.46 

827.65 

47.43 

.66 

1.00 

.40 

-7.44 

-7.16 

4.25 

9.12 

47.43 

4.47 

-.29 

-.07 

-.15 

50 

2.04 

-.72 

-.62 

6.07 

-2.69 

2.32 

.31 

-.16 

.05 

-.72 

.70 

-.91 

-2.69 

1.52 

-.44 

-.22 

.23 

-.01 

-.62 

-.91 

4.18 

2.32 

-.44 

3.03 

-.08 

.02 

-.08 

40 

_ 

.68 

1.32 

1.78 

.58 

-1.05 

-1.58 

.11 

.27 

-.06 

1.32 

4.55 

2.41 

-1.05 

12.35 

-1.59 

.20 

-.07 

-.13 

1.78 

2.41 

5.36 

-1.58 

-1.59 

6.87 

-.05 

.22 

.11 

50 

1 

_ 

23.19 

22.18 

16.56 

1 

9.27 

-.90 

.62 

.38 

.07 

-.02 

22.18 

33.81 

15.26 

I 

-.90 

7.03 

•1.82 

-.03 

.27 

-.01 

16.56 

15.26 

27.68 

1 

.62 

-1.82 

3.29 

-.32 

-.04 

-.00 

45 

r 

- 

2.93 

.38 

-.86 

| 

.83 

.73 

-1.86  I 

.19 

.07 

-.51 

.38 

.63 

-1.82 

I 

.73 

1.35 

-.69 

.03 

.27 

.00 

-.86 

-1.82 

7.55 

| 

-1.86 

-.69 

5.81  I 

.03 

-.03 

-.11 

j 

j 

. 

37 

- 

1.00 

-1.27 

1.11 

.99 

1.06 

-2.03 

.00 

.11 

.10 

-1.27 

2.91 

-.18 

1.  )6 

1.77 

-1.38 

.22 

.13 

-.35 

1.11 

-.18 

2.81 

-2.03 

-1.38 

5.54 

.04 

.14 

.06 

. 

H 


3.12 

-.06 

.09 

-.06 

9.75 

-2.85 

.09 

-2.85 

7.07 

2.49 

-.95 

1.48 

-.95 

5.86 

-1.81 

1.48 

-1.81 

7.06 

10.21 

5.64 

-4.12 

5.64 

5.19 

-1.34 

-4.12 

-1.34 

9.99 

5.92 

-2.69 

1.66 

-2.69 

1.75 

-1.36 

1.66 

-1.36 

7.74 

1.12 

-.75 

.42 

-.75 

17.94 

.67 

.42 

.67 

10.88 

21.40 

20.37 

20.41 

20.37 

32.48 

14.03 

20.41 

14.03 

31.01 

3.17 

1.00 

-1.55  ' 

1.00 

1 .48 

-2.47 

-1.55 

-2.47 

14.79 

- 

1.99 

-.69 

-1.19 

-.69 

4.10 

-1.36 

-1.19 

-1.36 

7.84 

.74 

1.71 

.93  : 

!  .69 

.05 

-2.72  ;  j 

.08 

-.08 

.01 

1.71 

8.30 

-2.18 

I  -05 

2.17 

-2-77  .  | 

-  .17 

-.06 

.40 

.93 

-2.18 

6.07  i 

i  -2.72 

-2.77 

14.07  i  ! 

.10 

-.04 

-.09 

1.31  2.16  -2.12 

2.16  10.97  -5 .03 

-2.12  -5.93  21.89 


33 


37 


!  1.94 

.24 

-.80 

i  !  .91 

-1.33 

.92  j 

-.09 

.66 

.09 

| 

3.10 

-1.37 

-1.18  i 

I  -24 

.50 

-1.60  1 

1  -1.33 

2.91 

-.31  j 

-.07 

.40 

-■01 

-1.37 

2.45 

-.30 

1  -.80 

-1.60 

7.16 

.92 

-.31 

2.46  1 

.49 

-1.39 

.25 

-1.18 

-.30 

7.49 

L 

r 

. 

L 

j  1 

1  | 

J 

- 

J 

[  2.17 

-3.28 

2.44 

1.56 

2.19 

-2.19  i 

-.21 

-.19 

.24 

1  4.52 

-1.28 

.11 

-3.28 

5.61 

-2.96 

2.19 

3.79 

-2.00  i 

.26 

.23 

-.32 

-1.28 

7.32 

-3.01 

1  2.44 

-2.96 

4.85  ; 

-2.19 

-2.00 

4.79  ; 

1  -.16 

-.11 

.19 

.11 

-3.01 

7.81  | 

J 

L 

J 

. 

. 

. 

J 

50 


0 


0 


.41 

.69 

.52 


-.69 

.52 

j  1.49 

2.06 

-2.61 

I  -.25 

-.07 

.19 

2.87 

1.15 

1  7 .  ‘.’6 

3.90 

-2.37 

.12 

.10 

-  .03 

1.15 

3.68 

j  -2.61 

-2.37 

6.33  1 

•  .27 

.05 

.34 

2.29  1.38  -2.10  ; 

1.38  6.12  -1.31  I 

-2.  0  -1.31  6.74  i 


37 


Date  T 

Ox 

Dy 

39 

• 

2.04  1.11  .17 

1.83  2.21 

1.11  1.14  1.39 

2.2!  3.08 

1 

.17  1.39  3.87 

-.60  -1.51 

55 

3.17  -1.58  2.08 

2.74  1.09 

-1.58  4.42  -2.51 

1.09  2.82 

2.08  -2.51  4.36 

-1.15  .09 

50  r 


54 


21  r 


46  r 


6/6/89  19 

MK-27 
C 


50 


33 


35 


35 


51 


1.04 

-1.43 

3.91 


41 


50  r 


-1.43  3.91 
12.39  -4.40 
-4.40  15.02 


1.80 

2.60 

2.80 


22.47 
!  2.34 

44.62 


.59 

-.42 

1.83 


2.14 

-2.12 

2.08 


1.92 

-1.31 

-.15 


814.12 

-28.74 

-10.84 


1.70 

-2.22 

2.57 


.54 

-.71 

1.51 


.50 

.02 

1.25 


.59  .76  .92 
.76  4.47  -2.17 
.92  -2.17  6.72 


8.84  25.18  -1.76 
25.18  107.44  .08 
-1.76  .08  4.53 


2.60  2.80 
5.00  2.22 

2.22  7.20 


2.34  44.62 

1.20  2.35 

2.35  108.62 


-.42  1.83 
4.12  -1.34 
•1.34  6.30 


-2.12  2.08 
3.03  -1.37 

-1.37  4.66 


-1.31  -.15 
1.57  -1.49 
-1.49  5.17 


-28.74  -10.84 
9.11  1.42 

1.42  1.63 


-2.22  2.57 
3.77  -1  75 
-1.75  7.55 


-.71  1.51 
3.71  -1.40 
-1.40  4.38 


.02  1.25 
2.50  -1.64 
-1.64  5.59 


-.60 

-1.51 

1.96 


-1.15 

.09 

4.54 


.45  -.05  -1.64 
-.05  1.96  -.06 
-1.64  -.06  6.26 


.56 

.47 

.76 


J  L 


.76 

-1.62 

-.26 


|  20.54 
|  2.78 

43.31 

r 

I  -27 

I  -.01 

-.71 


.69 

1.37 

.37 


1.50 

-.79 

.96 


2.78  43.31 
2.26  7.06 
7.06  108.13 


-.01  -.71 

2.32  .19 

.19  2. 


.79 

.76 

.04 


822.25 

-23.79 

1.99 


1.14 

-.06 

-4.60 


1 .04 
-.80 
-1.73 


-.06 

5.61 

1.32 


-4.60 

1.32 

18.88 


-.47  -.76 

2.03  .15 

.15  2.03 


10.11  28.52  .14 

28.52  111.27  -2.10 

.14  -2.10  1.66 


-1.62  -.26 
5.79  -2.58 
-2.58  4.46 


.96  j 
-.04  ! 
1.80 


23.79  1.99 

7.32  -.89 

-.89  4.39 


1.83  3.08  -1.64 
3.08  6.71  -1.07 
-1.64  -1.07  5.50 


.15 

.22 

.23 


.68 

.13 

-.01 


.03 

-.16 

.00 


.19 

.17 

-.31 


.75 

.93 

-.11 


-.01 

.02 

-.09 


1  f 


42 


1.37  .37] 
3.39  1.97  j 
1.97  2.81  ! 


.92 

j  .16 
1  -’2 

.34 

-.06 

.21 


.35 

-.12 

.19 


.12 

-.42 

.50 


1.00 

-.38 

-.19 


-.16 

.09 

-.12 


1  r 


-.80  -1.73  | 

2.09  .52  | 

.52  4.55 

1 


.08 

.10 

.05 


.39 

-.09 

.09 


Rxy 

.01 

.22 

.22 


.02 

.12 

.03 


.03 

.11 

.06 


.05 

.07 

.12 


.84 

.99 

.13 


.24 

.20 

.01 


.62 

.48 

.41 


.10 

.08 

.07 


.47 

.09 

.39 


.36 

.60 

.25 


.28 

.73 

.12 


.27 

.16 

.13 


.13 

.08 

.12 


.09 

.01 

.03 


.14 

.11 

.20 


.10 

.37 

.03 


.05 

.10 

.05 


.02 

.07 

.05 


.00 

.08 

.36 


.35 

.21 

.22 


.99 

.25 

.95 


1  [ 


J  L 


.14 

.13 

.13 


.07 

.16 

.09 


.15 

.07 

.31 


.09 

.14 

.08 


.01 

.00 

.07 


.18 

.15 

.12 


.03  1 
.03  j 
.14  | 

J 


3.37 

.91 

-4.02 


.59 

-.26 

1.05 


1.97 

-1.12 

2.09 


3.00 

-1.04 

.29 


.73 

.96 

-1.07 


4.12 

1.26 

1.25 


1.81 

-1.20 

-3.35 


.97 

-.54 

-.55 


.91 
1 .88 
.81 


-.26 

5.94 

-1.41 


-1.12 

6.99 

-1.20 


-1.04 

1.03 

-.98 


-1.20 

8.60 

-.40 


-.54 

4.53 

-1.10 


3.28 

2.7* 

-1.26 

2.75 

3.39 

-.81 

J 

-1.26 

-.81 

4.73 

1.93 

-.82 

1.35 

-.82 

6.37 

-3.56 

1.35 

-3.56 

9.16 

1.45 

-1.31 

2.40 

-1.31 

13.26 

•5.50 

2.40 

-5.50 

22.23 

.94 

.02 

.54 

.02 

6.92 

-2.15 

.54 

-2.15 

8.38 

1 

4.84 

-3.15 

-1.18 

| 

-3.15 

2.49 

1.54 

J 

-1.18 

1.54 

8.15 

2.58 

1.37 

1.99 

1.37 

12.90 

-1.62 

1.99 

-1.62 

9.19 

-4.02 

.81 

11.69 


1.05 

-1.41 

9.75 


2.09 

-1.20 

6.84 


.29 

-.98 

5.09 


.96  -1.07 
4.53  -.74 
-.74  6.45 


1.26  1.25 
8.85  -2.25 
-2.25  13.99 


-3.35 

-.40 

21.08 


-.55 

-1.10 

8.74 


38 


Date 

6/6/89 

MK-27 


C 


C 


cc 


7/21/88 

MK-27 


T 

50 


24 


41 


50 


50 


32 


27 


52 


28 


21 


50 


0 


17 


20 


Dx 

Dy 

Rxy 

M 

5.64 

1.24 

-.54 

1.87 

1.04 

-2.32 

.18 

-.29 

.05 

6.32 

2.39 

-2.57 

1.24 

.56 

-1.15 

1.04 

2.89 

2.46 

.05 

-.02 

-.04 

2.39 

3.50 

1.85 

-.54 

-1.15 

5.18 

-2.32 

2.46 

12.77 

-.05 

-.15 

.08 

-2.57 

1.85 

16.59 

. 

2.11 

-3.07 

1.33 

.96 

.77 

-2.12 

-.01 

-.16 

-.60 

3.11 

-2.01 

.03 

* 

-3.07 

4.83 

-1.40 

.77 

1.70 

.46 

.03 

.22 

.80 

-2.01 

5.29 

-2.59 

1.33 

-1.40 

2.33 

-2.12 

.46 

9.15 

.02 

.01 

-.09 

.03 

-2.59 

12.28 

.94 

-1.28 

1.88 

1.16 

.07 

-3.08 

.24 

.40 

.10 

1.60 

-2.39 

-1.00 

’ 

-1.28 

4.44 

-3.49 

.07 

5.25 

4.12 

.13 

-.12 

-.09 

-2.39 

10.89 

-.41 

1.88 

-3.49 

6.17 

-3.08 

. 

4.12 

14.24 

-.13 

.19 

.17 

-1.00 

-.41 

17.31 

. 

5.12 

-.37 

-3.60 

5.36 

-4.96 

3.52 

-.18 

.43 

.07 

12.33 

-6.79 

-.61 

■ 

-.37 

.92 

-1.39 

-4.96 

5.34 

-3.31 

.10 

-.03 

-.08 

-6.79 

6.39 

-4.00 

-3.60 

-1.39 

9.36 

3.52 

-3.31 

5.96 

. 

.01 

-.07 

-.11 

-.61 

-4.00 

16.90 

1.08 

.23 

2.28 

1.18 

1.03 

-1.97 

.56 

-.08 

-.08 

.99 

1.35 

.49 

1 

.23 

6.23 

1.96 

1.03 

8.48 

-5.15 

.04 

.27 

-.11 

1.35 

10.73 

-1.80 

2.28 

1.96 

9.27 

-1.97 

-5.15 

6.93 

.03 

-.06 

.05 

.49 

-1.80 

15.45 

J 

.75 

1.08 

.64' 

1.06 

-.70 

-1.98 

.14 

.19 

.12 

1.57 

.06 

-1.65 

1 

1.08 

3.68 

-1.80 

-.70 

2.14 

.89 

-.04 

.26 

.23 

.06 

4.37 

-1.68 

.64 

-1.80 

4.79 

-1.98 

.89 

5.98 

.02 

.01 

-.03 

-1.65 

-1.68 

11.09 

J 

50.78 

-108.04 

2.47 

59.74 

119.72 

4.42 

.96 

-.86 

.12 

5.24 

3.51 

2.24 

1 

108.04 

270.28 

.89 

119.72 

264.70 

-11.00 

1.06 

.99 

-.21 

3.51 

4.18 

3.47 

2.47 

.89 

9.99 

. 

4.42 

-11.00 

2.25 

.15 

-.12 

.11 

2.24 

3.47 

11.23 

5.46 

.27 

1.33 

r 

2.42 

1.17 

-.51 

.54 

.34 

-.03 

3.93 

.76 

-.88 

1 

.27 

.59 

-1.27 

1.17 

1.30 

.33 

.10 

.07 

-.10 

.76 

1  .T 

-1.12 

1.33 

-1.27 

11.53 

-.51 

.33 

8.84 

.27 

.25 

.51 

-.88 

-1.12 

10.15 

J 

1.17 

-1.54 

1.25 

.95 

.92 

-2.17 

.17 

.25 

-.10 

1.77 

-.79 

-1.15 

1 

-1.54 

2.42 

-.81 

.92 

1.66 

-.59 

-.15 

-.33 

-.11 

-.79 

5.42 

-1.31 

1.25 

-.81 

3.27 

-2.17 

-.59 

8.41 

.13 

.04 

-.35 

-1.15 

-1.31 

15.38 

. 

.65 

-1.23 

1.63 

.18 

.01 

-32 

.38 

.15 

-62 

.57 

-1.32 

.51 

-1.23 

6.29 

-4.26 

.01 

2.77 

.73 

-.04 

-.20 

-.16 

-1.32 

10.68 

-3.03 

1.63 

-4.26 

5.06 

-.32 

.73 

2.80 

.16 

.02 

.26 

.51 

-3.03 

5.88 

2.77 

-.50 

-.58 

2.26 

-1.81 

* 

1.83 

-.08 

.33 

-.13 

5.43 

-2.67 

.99 

-i 

-.50 

.36 

-.78 

-1.81 

1.82 

-1.27 

.03 

-.06 

-.01 

-2.67 

2.29 

-1.77 

-.58 

-.78 

4.50 

1.83 

-1.27 

2.28 

.27 

-.28 

.26 

.99 

-1.77 

5.13 

- 

.15 

-.27 

.73 

... 

.26 

-.05 

-.75 

.46 

-.10 

-.22 

.22 

.44 

-.36 

- 

-.27 

7.16 

-2.43 

-.05 

2.92 

-.24 

-.94 

.33 

.41 

.44 

7.03 

-3.21 

.73 

-2.43 

4.14 

-.75 

-.24 

2.28 

.99 

-.22 

-.51 

-.36 

-3.21 

9.57 

.87 

.20 

-.63 

[  '7i 

1.01 

-.65  ] 

.18 

.53 

-.27] 

1.32 

.73 

-.78 

1 

.20 

.78 

.95 

1  1.01 

2.61 

-1.97  1 

.05 

.58 

-.16  ! 

.73 

1.73 

-.62 

-.63 

.95 

3.95 

1  -.65 

L 

-1.97 

2.87  ' 

-  .02 

.07 

-  .12 

-.78 

-.62 

7.66 

39 


Date 


T 

50 


Ox 

4.19  -2.29  -2.04 
-2.29  1.66  2.19 
-2.04  2.19  4.80 


Dy 


Rxy 


2.66  -1.69  .90 

-1.69  1.49  -1.46 

.90  -1.46  2.85 


.34  -.23 
-.09  .06 
-.02  .01 


.00 

.09 

.14 


M 

4.61  -3.26  -1.08 
-3.26  2.97  .46 
-1.08  .46  6.60 


7/21/88  40 

MK-30 


117.59  -5.91  8.41 

-5.91  .67  2.50 

8.41  2.50  212.64 


114.89  -1.28  -12.84 
-1.28  3.71  -2.88 
-12.84  -2.88  233.43 


.98  -.34  -.01 
-.31  .17  .15 
-.02  -.44  .97 


5.55  2.32  -1.02 
2.32  3.85  .95 
-1.02  .95  14.68 


44.27 

-3.88 

5.53 

41.15 

4.05 

4.59 

.98 

.29 

.50 

-3.88 

3.57 

-.96 

4.05 

3.76 

2.17 

-.19 

-.36 

.14 

5.53 

-.96 

3.29 

4.59 

2.17 

10.53 

.22 

.24 

.22 

1.21 

-2.57 

1.78 

.27 

.53 

-.39 

.03 

-.04 

-.14 

-2.57 

6.27 

-3.16 

.53 

2.30 

.54 

-.01 

-.06 

-.03 

1.78 

-3.16 

3.51 

-.39 

.54 

2.18 

-.07 

-.18 

-.20 

1.51 

-.90 

-.35 

-.90 

9.96 

-.77 

-.35 

-.77 

11.23 

1 

.45 

-1.96 

1.64 

-1 

.96 

9.02 

-1.86 

1 

.64 

-1.86 

6.83 

2.29  -.21  -1.13 
-.21  .42  -1.06 
-1.13  -1.06  4.76 


1.99  -2.02  1.36 
-2.02  2.63  -.78 
1.36  -.78  2.46 


■ 

-.19 

.40 

-.07 

.06 

.04 

-.02 

.13 

-.57 

.17 

5.09  -2.73  .12 

-2.73  2.96  -1.19 

.12  -1.19  6.05 


119.29 

216.18 

5.06 

111.95 

216.93 

6.56 

216.18 

418.89 

1.31 

216.93 

439.23 

16.56 

5.06 

1.31 

6.37 

6.56 

16.56 

6.16 

99 

1.03 

.32 

92 

.99 

.36 

12 

.09 

-.03 

1.54  .44  -.19 

.44  9.76  -5.47 

-.19  -5.47  12.94 


.32 

.17 

.93 

.17 

6.96 

-3.21 

.93 

-3.21 

5.91 

1.76 

.77 

.34 

.77 

.90 

2.02 

.34 

2.02 

8.61 

.23 

.01 

-.62 

.01 

4.05 

-1.89 

-.62 

-1.89 

3.45 

.49 

1.05 

-.35 

1.05 

3.50 

-2.20 

-.35 

-2.20 

2.29 

.25 

-.16 

.03 

.28 

-.08 

-.20 

-.15 

-.17 

.14 

.11 

.03 

.07 

.05 

.07 

.00 

.01 

.05 

.06 

.41 

.06 

.43 

.06 

11.88 

-3.28 

.43 

-3.28 

8.08 

2. 

.04 

1.66 

-.18 

1. 

.66 

4.16 

-.28 

.18 

-.28 

10.36 

.88 

.97 

-1 

11 

.97 

5.98 

3 

29 

-1.11 

3.29 

7 

03 

.48 

-1.29 

1 

13 

-1.29 

6.55 

-1 

58 

1.13 

-1.58 

3 

87 

72 

.53 

.73 

53 

4.13 

-4.00 

73 

-4.00 

7.96 

1 

02 

-1 

68 

2.18 

-1 

68 

3 

92 

-1.97 

2 

18 

-1 

97 

7.12 

7.31 

-1 

41 

-4 

27 

-1.41 

1 

19 

1 

02 

-4.27 

1 

02 

5 

58 

1 

.77 

2.51 

2.26 

2 

.51 

4.16 

2.20 

2. 

.26 

2.20 

4.83 

5.41 

1.21 

-1.95 

.13 

.00 

.00 

5.71 

2.82 

-2.01 

1.21 

.58 

-1.61 

-.93 

-.16 

.59 

2.82 

7.15 

.86 

-1.95 

•1.61 

6.03 

-.47 

-.06 

.17 

-2.01 

.86 

10.89 

.90 

2.00 

-2.65 

.03 

-.03 

-.01 

1.34 

.45 

-1.65 

2.00 

7.22 

-4.99 

.18 

.05 

-.17 

.45 

13.08 

-5.32 

-2.65 

-4.99 

8.85 

.07 

-.04 

-1.65 

-5.32 

14.01 

.88 

-1.18 

-1.05 

.18 

-.30 

.34 

1.30 

.57 

-1.80 

-1.18 

9.74 

-11.81 

-.25 

.00 

.27 

.57 

13.83 

-17.34 

-1.05 

-11.81 

23.56 

.14 

-.08 

-.04 

-1.80 

-17.34 

32.58 

3.70 

5.55 

-5.68 

-.08 

.09 

.16 

5.03 

3.31 

-3.86 

5.55 

11.56 

-4.41 

.06 

.02 

-.05 

3.31 

15.25 

-7.42 

-5.68 

-4.41 

14.07 

-.13 

.20 

.30 

-3.86 

-7.42 

15.10 

3.59 

-2.61 

1.15 

.03 

-.03 

.11 

10.64 

-3.53 

-3.19 

-2.61 

2.49 

-.79 

-.10 

-.13 

.08 

-3.53 

4.15 

-.21 

1.15 

-.79 

2.42 

-.10 

.08 

.21 

-3.19 

-.21 

6.50 

1.31 

-3.38 

1.14 

.25 

-.54 

.77] 

2.31 

-  .16 

2.36 

-3.38 

9.70 

-4.71 

.14 

-.29 

.15 

-.16 

17.52 

-2.20 

1.14 

-4.71 

4.80 

.12 

-.29 

.16  ; 

2.36 

-2.20 

8.09 

j  l  J 
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APPENDIX  D.  AUTOCORRELATIONS 


Examples  of  autocorrelation  functions  of  residuals  R(h)/R(0),  (see  eq.  (41)) 
are  presented  in  Figure  D1  through  D6.  The  first  two  are  for  the  displacement 
process,  eq.  (40).  The  latter  four  are  for  the  three  components  of  displacement. 
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Array  4 


Array  5 


Figure  D3.  Autocorrelation  of  noise  components 
5/10/89  5  =  1.24 
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Figure  D4.  Autocorrelation  of  noise  components 
5/10/89  6  =  2.85 
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Downrange 


Downrange 


Array  1 


Array  2 


Downrange 


Downrange 


APPENDIX  E.  COMPUTER  SOURCE  CODES 


The  major  programs  needed  to  produce  Table  1  are  documented  here. 
Basically  there  are  three  steps: 

(i)  Extract  pertinent  crossover  data  from  T-files 

(ii)  Produce  the  covariance  matrices  of  residuals;  Dx,  Dy,  Dxy  and  M.  See 
equations  (38)  and  (39) 

(iii)  Develop  the  estimates,  eq.  (9)  and  (10).  Develop  supporting 
statistics,  eq  (17),  (18),  (19) 

Also  the  user  needs  access  to  Table  E.l,  the  coordinates  of  the  various 
arrays  at  the  Nanoose  range. 

Step  (i)  is  managed  by  the  program  KEYGATE,  written  in  FORTRAN  77, 
Miscosoft  optimizing  compiler  4.01.  This  program  reads  the  NUWES  T-files 
and  performs  a  series  of  gating  operations  in  order  to  produce  crossover  data 
in  the  required  seven  column  format.  It  will  prompt  the  user  for 

(a)  name  of  the  range  (e.g.,  Nanoose) 

(b)  Number  of  records  to  ignore  (e.g.,  31  for  bypassing  the  DVT 
information) 

(c)  The  target  vehicle  mode  (e  g.  7) 

First,  it  will  strip  out  all  data  other  than  mode  2  or  mode  7.  Second,  it 
identifies  all  point  counts  for  which  a  position  vector  is  available  for  two  or 
more  arrays.  Third,  it  reads  an  array  location  file  and  removes  data  that 
cannot  be  reasonably  identified  with  an  overlap  region.  Next,  it  organizes  the 
data  by  pairs  of  arrays.  Finally,  it  prompts  the  user  for  his  array  pair  and 
output  filename;  it  selects  and  records  all  the  data  of  the  desired  type. 
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Step  (ii)  is  contained  in  the  program  PRINCOM3.M.  This  is  a 
PCMATLAB  program  that  performs  the  eigen  analysis,  eq.  (32);  develops  the 
smoothed  fit,  eq.  (37);  and  the  covariance  matrices  (38)  and  (39).  The  latter  is 
placed  on  an  output  file  to  be  used  in  Step  (iii).  It  also  develops  the 
autocorrelation  sequence,  e.g.  (41)  of  the  displacement  process. 

A  slight  modification  of  this  program  can  be  made  to  develop 
autocorrelations  for  the  three  components  of  residuals  instead  of  for  the 
displacement  process. 

Step  (iii)  is  performed  by  the  program  TIMECOR,  a  FORTRAN  77  code. 
The  Microsoft  FORTRAN  optimizing  compiler  4.01  is  used.  The  user  should 
note  that  requirements  of  the  three  input  files: 

VELOCITY.DAT  is  a  two-column  data  set  containing  the  water  layer 
boundaries  and  the  sound  speeds  in  25  ft  increments 

TRPDOTRX.DAT  is  the  seven-column  crossover  data  output  of  KEYGATE 
to  which  has  been  appended  the  locations  of  the  two  arrays  as  the  ’ast 
record. 

MMATRIX.M  is  the  covariance  matrix,  M,  produced  by  PRINCOM.M  in 
Step  (ii). 

Since  the  ray  tracing  exit  layer  elevation  angles  are  not  contained  in 
T-files  it  is  necessary  to  reconstruct  them  in  order  to  compute  Equation  (5). 
This  is  accomplished  by  the  subroutine  TGEN.  The  methodology  of  TGEN  is 
explained  in  [3]  under  the  heading  of  ray  fitting. 

There  are  two  output  files  which,  taken  together,  contain  all  of  the 
information  indicated  in  Table  1. 


TABLE  El.  COORDINATES  OF  THE  NANOOSE  ARRAYS 


Array  Number 

Date  of  Survey 

Downrange 

Crossrange 

Depth 

0 

6/20/85 

12188.01 

-131.52 

-1295.33 

1 

6/20/85 

19463.16 

-174.99 

-1308.76 

2 

7/12/85 

26991.39 

-109.83 

-1323.25 

3 

1/07/88 

34505.10 

-80.76 

-1323.32 

4 

10/24/88 

42005.19 

-55.17 

-1318.28 

5 

6/20/85 

49497.00 

-25.23 

-1315.58 

6 

6/20/85 

56972.28 

-21.21 

-1308.50 

7 

7/30/85 

64680.66 

15.33 

-1353.39 

8 

11/16/88 

71969.73 

-29.28 

-1300.89 

9 

5/07/84 

3.00 

3.00 

1.00 

10 

3/12/84 

47100.00 

-3600.00 

-1300.00 

11 

7/18/85 

23173.89 

-6488.40 

-1312.09 

12 

6/20/85 

30731.25 

-6553.05 

-1312.90 

13 

6/20/85 

38213.61 

-6640.77 

-1323.05 

14 

6/20/85 

45647.07 

-6513.18 

-1324.78 

15 

6/19/85 

53249.43 

-6354.60 

-1316.66 

16 

9/13/85 

60859.74 

-6356.07 

-1313.42 

17 

6/16/87 

68217.93 

-6524.10 

-1313.43 

54 

2/02/88 

38029.95 

5401.98 

-1212.69 

55 

6/20/85 

45645.75 

6369.66 

-1188.12 

56 

7/30/85 

53180.13 

6417.96 

-1218.84 

57 

7/30/85 

60745.71 

6419.40 

-1088.24 

23 

6/20/85 

41605.14 

-12150.18 

-1268.23 

24 

4/17/89 

49572.00 

-12966.00 

-1300.00 

25 

10/24/88 

56993.79 

-12999.3? 

-1205.48 

26 

8/08/88 

64442.94 

-12971.04 

-1255.35 

27 

7/15/80 

22119.60 

-15908.70 

83.00 

28 

5/04/83 

45000.00 

1500.00 

-1350.00 

29 

2/02/79 

.00 

.00 

.00 

KEYGATE 


PROGRAM  KEYGATE 


Program  to  read  in  raw  data  from  Keyport  hydrophone  arrays,  segregate  it  by  mode,  and  throw  out 
unusable  records.  The  output  of  this  program  is  to  be  read  in  by  the  program  KEYMA1N. 

Modified  by  Colin  R.  Cooper  11/15/89 


INTECER'4  CRT,  KBD 
CHARACTER'25  DSNAME,  SITNAM 
CHARACTERS  TEMPI,  TEMP2,  TEMP3,  TEMP4 

PARAMETER  (KBD=5,  CRT  =6) 


WRITE(CRT,*)  1  Please  enter  the  name  of  your  input  file:  ' 
READ(K8D/(A)')  DSNAME 

WRITEfCRT,*)  '  Please  enter  the  name  of  the  range  configuration  file:', 
READ(KBD,'(A)')  SITNAM 

TEMPI  =  TEMPI  TMP' 

TEMP2  =  TEMP2.TMP' 

TEMP3  =  TEMP3.TMP' 

TEMP4  =  TEMP4.TMP' 

CALL  STRMOEKCRT,KBD,TEMPl, DSNAME) 

CALL  PAIR(CRT,KBD,TEMP1,TEMP2) 

CALL  RANGE(CRT,KBD,SITNAM,TEMP2,TEMP3) 

CALL  PAlR2(CRT,kBD,TEMP3,TEMP4) 

CALL  REC(CRT,KBD.TEMP4) 

WRITE(CRT,*)  Operation  complete.  KEYGATE  terminating... 

STOP 

END 


SUBROUTINE  STRMOD(CRT,KBD, TEMPI  ,DSN AME) _ 

Program  to  strip  all  modes  except  2's  and  7's  from  the  Keypon  data.  (2  indicates  target  ship,  7 
indicates  torpedo). 


CHARACTER  DO*2,  DSN'AME*25,  TEMP1S 
INTEGER  PC,  ARRAY,  CRT,  NHEAD 

OPEN(l,FILE=DSNAME,STATUb='OLD) 

VVRITEiCRT,')'  Ho"  manv  records  of  header  do  vou  want  to  strip  off  the  file'1', 
READ(KBD,*)  NHEAD 
DOll  1  =  1, NHEAD 
READfi.*) 

11  CONTINUE 
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WRITEfCRT,*)'  Input  mode  to  be  kept?' 

READfKBD,*)  NUM 

OPEN(2,FILE=TEMPl  ,ST  ATUS='NE  W') 

10  READ(1,100,END=50,ERR=40)PC,DO,X,Y,Z,ARRAY,MODE 
lF(DO.NE.'  ')GOTO  20 
IF(MODE  .NE.  NUM)  GOTO  20 
WRITE(2,1 1 0)PC,X,Y,Z,  ARRA  Y.MODE 
20  CONTINUE 
GOTO  10 

40  WRITE(CRT/*)'  THERE  WAS  AN  ERROR  IN  THE  FILE' 
50  CONTINUE 

100  FORMAT(I5/A2,1X,F7.1,2X,F7.1,2X/F7.1,30X,I2,2X,I1) 

110  FORMAT(1X,15,2X,3F10.1,2X,12,2X,I2) 

close  (UNrr=i) 

CLOSE(UNIT =2) 

RETURN 

END 


SUBROUTINE  PAIR(CRT,KBD,TEM PI  ,TEMP2) 


Program  to  pair  point  counst  after  the  data  has  been  gated  by  STRMOD.  Second  pass. 


DIMENSION  X(200),  Y(200),  Z(200) 

INTEGER‘4  PC(200),  ARRAY{200),  MODE(200),  CRT.  HOLD 
CHARACTER-9  TEMPI,  TEMP2 

OPEN(l,FILE=TEMPl,STATUS='OLD) 

OPEN(2,FILE=TEMP2,STATUS='NEW'l 


HOLD  =  0 
IREC  =  0 
NREC  =  0 
1  =  1 


Read  records  by  two's  to  compare  point  counts. 


READ(1,*,END=40,ERR=30)  PC(I),  X(I),  Y(I),  Z(I), 

+  ARRAY(I),  MODE(I) 

10  READ(1,*,END=40,ERR=30)  PC(I+1),  X(M),  Y<I+1),  Z0+1), 
+  ARRAY0+1),  MODEO+l) 

NREC  =  NREC  +  1 


IF(PC(I+1)  .EQ  HOLD)  THEN 

WRTE(2.100)  PC(I+1),  X(I-l),  Y(I+1),  Z(I+1), 
ARRAY(I-l),  MODE(l+l) 

HOLD  =  PC0+1) 

GO  TO  20 
END  IF 

IFfPC(I)  .EQ.  PCfl+l))  THEN 

WRTEC2.100)  PC(I),  X(I),  Y(I),  Z(I),  ARRAY(I),  MODE(I) 
WRTE(2,100)  PC(I+1),  X(I+1),  Y(I+1),  Z(M), 


i 


OPEN(l,FILE=TEMP2,STATUS='OLD') 
OPEN(2,FILE=SITNAM/STATUS='OLD') 
OPEN(3,FlLE=TEMP3,ST  ATUS='NEW') 


Read  site  configuration  into  CONFIG  array: 


NREC  =  0 
I  =  1 

10  READ(2,*,END=40,ERR=30)  CONFlG(l,l ),  CONFlG(I,2),  CONFIC(I,3), 
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+  CONFlC(I,4) 

NREC  =  NREC  +  1 
1  =  1  +  1 
GOTO  10 

30  WRITEfCRT,*)'  There  was  an  error  reading  the  config  file  in  record',NRB 
40  CONTINUE 
CLOSE(UNTT =2) 

NDREC  =  0 


Read  X,  Y,  Z,  and  ARRAY  from  input  data  file: 


45  RE  ADO  ,*,ERR=70,END=80)  PC,  X,  Y,  Z,  ARRAY,  MODE 
DO  501  =  1,NRE 
LX  =  CONFlG(l,l) 

LY  =  CONFICfU) 

LZ  =  CONFlG(I,3) 

CARRAY  =  INT(CONFlG(l,4)) 

MAXVAL  =  4700. 


Match  array  number  in  data  file  with  that  in  config.  file.  If  they  are  equal  compute  slant  range 
distance  (SRD1ST): 


IF( ARRAY  .EQ.  CARRAY)  THEN 

SRD1ST  =  SQRT(  (X  -  LX)~2  +  (Y  -  LY)"2  +  (Z  -  LZ)**2  ) 
1F(MAXVAL  .GE.  SRDIST)  THEN 

WRITE(3,100)  PC,  X,  Y,  Z,  ARRAY,  MODE,  SRDIST 
NDREC  =  NDREC  +  1 
END  IF 
END  IF 

50  CONTINUE 
GOTO  45 

70  WRITEfCRT,*)-  There  is  an  error  in  the  data  file.  ' 

80  CONTINUE 

CLOSE  (UNIT=1) 

CLOSE(UNIT=3) 

100  FORMAT  (3X,I5,3(3X,F1 0 . 1  ),2(3X,I2),3X,F8.2) 

RETURN 

END 


SUBROUTlNEPAlR2(CRT,KBD,TEMP3)TEMP4) 

Program  to  pair  point  counts  after  the  data  has  been  tested  by  RANGE.  Fourth  pass. 


DIMENSION  X(200),  Y(200),  Z(200),  SRDIST(200) 
INTEGERS  PC(200),  ARRAY(200),  MODE(200),  CRT,  HOLD 
CHARACTERS  TEMP3,  TEMP4 
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OPENd  ,FILE=TEMP3,ST  ATUS='OLD) 
OPEN(2,FILE=TEMP4,STATUS='NEW') 


HOLD  =  0 
IREC  =  0 
NREC  =  0 
1=1 


Read  records  by  two's  to  compare  point  counts. 


READd/*,END=40,ERR=30)  PCU),  Xd),  Y(l),  Z(I), 

+  ARRAYd),  MODEd),  SRDIST(I) 

10  READd,*,END=40,ERR=30)  PCd+1),  Xd+1),  Yd+1),  Zd+1), 

+  ARRAYd+1),  MODEd+1),  SRDISTd+1) 

NREC  =  NREC  +  1 

IF(PCd+l)  .EQ.  HOLD)  THEN 

WRITE(2,100)  PCd+1),  Xd+1),  Yd+1),  Zd+1), 

+  ARRAYd+1),  MODEd+1),  SRDISTd+1) 

HOLD  =  PCd+1) 

CO  TO  20 
END  IF 

IF(PC(I)  .EQ.  PCd+1))  THEN 

WRITE(2,100)  PC(I),  X(l),  Y(I),  Z(I),  ARRAY(I),  MODEd), 

+  RDISTd) 

WRITE(2,100)  PC(I*1),  Xd+1),  Y(!+l),  Z(I+1), 

+  ARRAYd+1),  MODEd+1),  SFD1ST(1+1) 

HOLD  =  PCd+1) 

IREC  =  IREC  +  1 
CO  TO  20 
END  IF 

IF(PC(I)  .NE.  PCd+1))  THEN 
PCd)  =  PCd+1) 

xo)  =  xd+i) 

Y(l)  =  Yd+1) 

Z(l)  =  Zd+l) 

ARRAY(l)  =  ARRAY(1+1) 

MODEd)  =  MODEd+1) 

SRDISTd)  =  SRDIST(1+1) 

1  =  1 
END  IF 

GOTO  10 
20  1  =  1+1 
GOTO  10 

30  WRITECCRT,*)'  THERE  IS  AN  ERROR  IN  THE  DATA  FILE  IN  REC.',NRE 
WRITE(CRT,*)  ...  OPERATION  TERMINATING  DUE  TO  ERROR.  ' 

STOP 

40  CONTINUE 

CL05E(UNTT=1) 

CLOSE(UNTT=2) 

100  FORMATdX,I5,2X,F7.1,2X,2F9.1,2X,12,2X,I2,2X,F8.2) 
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RETURN 

END 


SUBROUTINE  REC(CRT,KBD,TEMP4) 


Program  to  produce  the  final  gating  of  Keyport  hydrophone  array  test  data. 


INTEGER  PC,  ARRAY,  PCI,  PC2,  Al,  A2,  ARRAY1,  ARRAY2,  CRT 
INTEGER  PCA(IO),  ARRAYA(IO),  HOLD 
DIMENSION  XA(10),  YA(10),  ZA(10),  SRDISTA(IO) 

CHARACTER  OUTFlL*25,  ANS’l,  TEMP4*9,  TEMPI *9,RANGE*7 

RANGE='NANOOSE' 


99  WRITE(CRT,*)'  What  is  the  name  you  wish  to  give  to 
+  '  your  output  file?  ' 

READ(CRT,'(A)')  OUTFIL 

WRITE(CRT,’)'  Enter  the  numbers  of  the  arrays  to  be  paired:  ' 
REAEXKBD,*)  Al,  A2 

OPEN(l,FILE=TEMP4,STATUS='OLD) 

OPEN(2,FILE='TEMPl.DAT',STATUS='OLD') 

5  RE AD(1  ,*,END=8,ERR=8)  PC,  X,  Y,  Z,  ARRAY,  MODE,  SRDIST 
IF(( ARRAY  .EQ.  Al).OR.( ARRAY  .EQ.  A2))  THEN 
WRITE(2,100)  PC,  X,  Y,  Z,  ARRAY,SRD1ST 
END  IF 
GOTO  5 

8  CONTINUE 
CLOSE  (UNIT=1) 


Routine  to  pair  data  again: 


REWIND  2 

OPEN(4,FILE=TEMP2.DAT',STATUS=OLD) 


1  =  1 

I  FLAG  =  0 
FIRST  =  1 
HOLD  =  0 
M  =0 


9  REACX2,100,END=40)  PCA(I),  XAO),  YA(I),  ZA(I), 
+  ARRAYA(I),  SRDIST  AO) 

IF(FIRST  .EQ  1)  THEN 
RRST  =  0 
HOLD  =  PCA(I) 

END  IF 

IF(PCAO)  .EQ.  HOLD)  THEN 
1  =  1  +  1 
M  =  M  +  1 
GOTO  9 


ELSE 


In  cases  where  there  are  three  or  more  reports  for  a  given  point  count,  segregate  by  comparing 
SRDIST.  if  there  are  more  than  3  reports,  discard  all. 

IF  (M  .EQ.  3)  THEN 

IFCARRAYA(M)  .EQ.  ARRAY  A(M-l))  THEN 
IF  (ARRAY A(M)  .EQ.  ARRAYA(M-2))  THEN 
GOTO  50 
ELSE 

IF  (ABS(SRDISTA(M-2)-SRDISTA(M)>  .LT. 

+  ABS(SRDISTA(M-2)-SRDISTA(M-l)))  THEN 

WRITE(4,100)  PCA(M),  XA(M),  YA(M),  ZA(M), 

+  ARRAYA(M),  SRDISTA(M) 

WRITE(4,100)  PCA(M-2),  XA(M-2),  YA(M-2),  ZA(M-2), 

+  ARRAYA(M-2),  SRDISTA(M-2) 

IFLAG  =  1 
ELSE 

WRITE(4,100)  PCA(M-l),  XA(M-l),  YA(M-l),  ZA(M-l), 

+  ARRAYA(M-l),  SRDISTA(M-l) 

WRITE(4,100)  PCA(M-2),  XA(M-2),  YA(M-2),  ZA(M-2), 

+  ARRAYA(M-2),  SRDISTA(M-2) 

IFLAG  =  1 
END  IF 
END  IF 
ELSE 

IF  (ARRAY ACM)  .EQ.  ARRAY A(M-2)>  THEN 
IF  (ABS(SRDlSTA(M-l)-SRDISTA(M))  .LT. 

+  ABS(SRDISTA(M-l)-SRDlSTA(M-2)))  THEN 

WRITE(4,100)  PCA(M),  XA(M),  YA(M),  ZA(M), 

+  ARRAYA(M),  SRDISTA(M) 

WRITE(4.100)  PCA(M-l),  XA(M-l),  YA(M-l),  ZA(M-l), 

+  ARRAYA(M-l),  SRDISTA(M-l) 

IFLAG  =  1 
ELSE 

WRITE(4,100)  PCA(M-l),  XA(M-l),  YA(M-l),  ZA(M-l), 

+  ARRAYA(M-l),  SRDISTACM-D 

WRITE(4,100)  PCACM-2),  XA(M-2),  YA(M-2),  ZA(M-2), 

+  ARRAYACM-2),  SRDISTA(M-2) 

IFLAG  =  1 
END  IF 
ELSE 

IF  (ABS(SRDISTA(M)-SRDISTA(M-l))  .LT. 

+  ABS(SRDlSTA(M)-SRDISTA(M-2)))  THEN 

WRITE(4,100)  PCA(M),  XA(M),  YA(M),  ZA(M), 

+  ARRAYA(M),  SRDISTA(M) 

WRITE(4,100)  PCA(M-l),  XA(M-l),  YA(M-l),  ZA(M-l), 

+  ARRAYA(M-l),  SRDlSTA(M-l) 

IFLAG  =  1 
ELSE 

VVRITE(4,100)  PCA(M),  XA(M>,  YA(M),  ZA(M), 

+  ARRAYA(M),  SRDIST A(M) 

WRITE(4,100)  PCA(M-2),  XACM-2),  YA(M-2),  ZA(M-2), 

+  ARRAYA(M-2),  SRDISTACM-2) 

IFLAG  =  1 
END  IF 
END  IF 
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END  IF 


ELSE 

IF  ((M  .EQ.  2)  .AND.  (ARRAYA(M)  .NE.  ARRAYACM-l)))  THEN 
WRITE(4,100)  PCA(M),  XA(M),  YA(M),  ZA(M),ARRAYA(M), 
+  SRDISTA(M) 

WRITE(4,100)  PCACM-l),  XACM-l),  YA(M-l),  ZA(M-l), 

+  ARRAYACM-l),  SRDISTACM-2) 

END  IF 
END  IF 

50  CONTINUE 

25  CONTINUE 

IFLAG  =  0 
PCA(l)  =  PCA(I) 

XA(1)  =  XA(I) 

YA(1)  =  YA(I) 

ZA(1)  =  ZA(I) 

ARRAYA(l)  =  ARRAYA(I) 

SRDISTA(l)  =  SRDISTA(I) 

1  =  2 
M  =  1 

HOLD  =  PCA(l) 

D045L=  2,10 
PCA(L)  =  0 
XA(L)  =  0 
YA(L)  =  0 
ZA(L)  =  0 
ARRAY  A(L)  =  0 
SRDISTA(L)  =  0 
45  CONTINUE 

END  IF 

GOTO  9 

40  CONTINUE 


NREC  =  0 


CLOSE  (UNIT=4) 

OPEN(7,FILE=TEMP2.DAT',STATUS='OLD) 

OPEN(9,FlLE=OUTFIL,STATUS=NEW) 

WRITE(9,300)  RANGE,  Al,  A2 


Read  in  array  data  in  two  record  pairs: 


10  READ(7,100,END=60,ERR=70)  PC3,  X),  Yl,  Zl,  ARRAY),  SRDJST1 
READ(7,100,END=60,ERR=70)  PC2,  X2,  Y2,  Z2,  ARRAY2,  SRD1ST2 

IF(ARRAY)  .EQ.  Al  .AND.  ARRAY2  EQ.  A2)  THEN 


r 


1/  arrays  are  in  specified  order  (e.g.  4,5): 


WRrTE{9,200)  PCI,  XI,  Yl,  Zl,  X2,  Y2,  Z2 
END  IF 

IF(ARRAY1  .EQ.  A2  .AND.  ARRAY2  .EQ.  Al)  THEN 


If  arrays  are  in  reverse  order  (eg.  5,4): 


WRITE(9,200)  PCI,  X2,  Y2,  Z2,  XI,  Yl,  Zl 
END  IF 


Increment  record  counter: 


NREC  =  NREC  +  1 
GOTO  10 


70  WRITEfCRT,*)'  There  is  a  bad  record  in  the  file.’ 
60  CONTINUE 

CLOSE  (UNIT=2) 

CLOSEfUNIT =7) 

CLOSE(UNIT=9) 

100  FORMAT(I5,3(2X,F8.1),2X,I2,2X,F8.2) 

200  FORMAT(2X,15,lX,6(2X,Fll.l)) 

300  FORMAT(16X,A10,2X,I2,3X,I2) 


WRITEtCRT,*)  '  Do  you  want  to  try  another  array  pair?  (Y/N)' 

READ(KBD,’<AV)  ANS 

IF(  ANS  .EQ.  V  .OR.  ANS  .EQ.  y  )  GO  TO  99 

RETURN 

END 


function!  Autocorx,Autocoryl^princom3(fname) 

3/13/90 

PRINCOM3.M  will  evaluate  the  principal  components  of  the  passed  data  file.  Computes  the 
autocorrelation  of  the  displacement  process. 


iname  =  ('d:\mdataV  fname  ’.out']; 

evaUJ'load  '  iname  ]), 
eval(['keyout  =  '  fname 
cval([  clear  '  fname]); 
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X=keyout(:/5:7); 

Y=keyout(:,2:4); 

(len,n]=size(keyout); 


Compute  averages  and  covariance  matices  for  the  tracks  produced  by  each  of  the  arrays. 

Xav=onesflen,l)*sum(X)/len  ; 

Yav=onesflen,l)*sum(Y)/len  ; 

XX=X-Xav ; 

YY=Y-Yav ; 

ipc=keyout(:,l); 

CX=cov(XX); 

CY=cov(YY); 


Develop  the  eigen  analysis.  Choose  largest  eigen  value. 

[WX,DX)=eig(CX)  ; 

[  WY,DY]  =eig(CY)  ; 

ifDXO,l)>DX(2,2) 

if  DX(1,1)>DX(3,3) 
bc-1; 
end 

elseif  DX(2,2)>DX(3,3) 
lx  =2  ; 

end 

if  DX(3,3)>DX(1,1) 

if  DX(3,3)>DX(2,2) 
lx=3 ; 
end 

end 

ifDY(l,l)>DY(2,2) 

if  DY(1/1)>DY(3,3) 
ly=i; 
end 

elseif  DY(2,2)>DY(3,3) 
ly=2  ; 

end 

if  DY(3,3)>DY(1,1) 

if  DY(3r3)>DY(2,2) 

Iy=3; 

end 

end 

PX=XX'WX; 

PY = YY*  WY ; 

ux=PX(:,lx); 

uy=PY(:,ly); 


Modify  data  projection  onto  the  first  principal  component  to  account  for  constant  speed. 


ipcav=sum(ipc)/length(ipc); 
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mx=(sum(ux.*(ipc-ipcav)))/(sum((ipc-ipcav).*(ipc-ipcav))) ; 
my=(sum(uy.*(ipc-ipcav)))/(sum((ipc-ipcav).*(ipc-ipcav))) ; 

ax=(sum(ux)/length(ux))  -  (mx’ipcav); 
ay=(sum(uy)/length(uy))  -  (my*ipcav); 

uhx=(mx*ipc)+ax; 

uhy=(my*ipc)+ay; 

PPX=zeros(len,3); 

PPY=zeros(len,3); 

PPX(:,lx)=uhx; 

PPY(:,ly)=uhy; 


Compute  the  straight  line  tracks  in  the  original  coordinate  system. 


XSL=(PPX*WX')+Xav; 

YSL=(PPY*WY')+Yav; 


Develop  the  displacement  process  of  residuals. 


distx  =  sqit((X(:,l)  -  XSL(:,1)).A2  ♦  (X(:,2)  -  XSL(:,2)).A2  +  (X(:,3)  -  XSL(:,3)).*2); 
disty  =  sqrt((Y(:,l)  -  YSL(:,1)).A2  +  (Y(:,2)  -  YSL(:,2)).A2  +  (Y(:,3)  -  YSLf:,3)).A2); 

davx  =  sum(distx)/length(distx); 
davy  =  sum(disty)/length(disty); 

vx  =  distx  -  davx; 
vy  =  disty  -  davy; 


Compute  the  auto  correlations. 


Autocorx  =  zeros(length(vxl,l); 

Autocory  =  zeros(length(vy),l); 

for  k  =  l:length(vy) 

for  i  =  l:length(vv)  -  k  +  1 

Autocorx(k)  =  Autocorxfk)  +  vx(i)*vx(i  +  k-1); 

Autocory(k)  =  Autocoryfk)  +  vy(i)*vy(i  +  k-1); 
end 

Autocorx(k)  =  Autocorx(k)/sum(vx./'2); 

Autocoryfk)  =  Autocory(k)/sum(vy,A2); 

end 

dg 

subplot(211) 

plot(Y(:,l),Y(:,2),'or',YSL(:,l),YSL(:,2),'-r',X(:,l),X(:,2),'xg',XSL(:,l),XSL(:,2),'-g') 
title([fname  '  -  Corrected  Tracks  with  Principal  Components'!) 
xlabclCDown  Range  ), ylabelCCross  Range’) 
subplot(212) 

plot(Y(:,l),Y(:,3),'or’,YSL(:,l),YSL(:,3),'-r',X(:,l),X(:,3),'xg’,  XSL(.,l),XSL(:,3),'-g') 
titlef'Corrected  Tracks  with  Principal  Components') 
xlabclCDown  Range'),ylabel( 'Depth') 
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DIMENSION  IPCd  50) 

CHARACTER*30  LINE 

REA  L*8  Y1  (1 50.2). Y2(l  50,2), Y3(l  50,2),' YCl  (1 50,2), YC2(  1 50,2) 
REAL*8  LL(55),G(55),W (55), A2,P1  ,P2,V0,Vl ,TT1ME(1 50,2) 
REAL*8  THETA(150,2),DEPTH(55),PHia50,2),V(150,2),DZ 
REAL*8  Bl  (1 50,2), B2(150, 2),B3(1 50,2), YC3(1 50,2) 

REAL‘8  DEN,GG,SSB,SSR,TD,TM,SIG,DD,TEMP 
REA  1*8  MEQ,CEQ(150),DELTA,T1,TNOT,D1,D2,MSR,MSB 
REAL’S  DUMl,DUM2,DUM3,DUM4,DUM5,DUM6,DUM7 
REAL*8  M(3,3),C(150,3),VDEL,SDEL,VM,SDVM.COVDM 
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REAL*8  DELTAN,MEQN,R,CX(3,3),CY(3,3),CXY(3,3),RXY(3,3),T 

INTEGERS  TEE 


Read  in  the  data  from  the  data  files. 


OPEN(UNlT=2,FILE='VELOCITY.DAT',STATUS='OLD') 

OPENCUNn^FILE^RPDOTRX.DAT^TATUS^OLD) 

OPEN(UNnT=10/FILE=’MMATRIX.M,,STATUS='OLD) 

OPEN(UNIT=ll,FILE='OUTPUTl.COV’,STATUS='OLD) 

OPEN(UNIT=12,FILE='OUTPUT2.COV',STATUS='OLD‘) 

RE  AD(10,1 20)CX(1 ,1),CX(2,1  ),CX(3,1  ),CX(1 ,2),CX(2.2), 

*  CX(3,2),CX(1,3),CX(2,3),CX(3,3) 
REAEX10,120)CY(l,l)/CY(il),CY(3,l)/CY(U)(CY(2^), 

•  CY(3,2),CY(1,3),CY(2,3),CY(3,3) 

RE  AD(10,1 20)CXY(1 ,1  ),CXY(2,1),CXY(3,1  ),CXY(1 ,2),CXY(2,2), 

*  CXY(3,2),CXY(1,3),CXY(2,3),CXY(3,3) 

READG0,124)T 

TEE  =  IDINTCT) 

IF  (TEE.EQ.O)  THEN 
IX)  2  1  =  1,3 

RXY(1,I)  =  O.ODO 
RXY(2,1)  =  O.ODO 
RXY(3,I)  =  O.ODO 

2  CONTINUE 
ELSE 

DO  4  I  =  1,3 
D03  J  =  1,3 

RXY(J,I)=CXY(J,I)/DSQRT(CX(I,I)*CY(J,J)) 

M(J,D  =  CX(],1)  +  CY(J.l)  -  CXY(IJ)  -  CXY(j,I) 

3  CONTINUE 

4  CONTINUE 
ENDIF 

5  READ(12,'(A)',END=7)LINE 
GOTO  5 

7  BACKSPACE  12 
WRITE(12,122)TEE 
D08I  =  13 

8  WRITE(12.121K:X(I,1),CX(I,2),CX(1,3),CY(I,1),CY(I,2), 

•  CYd,3),RXYd,l),RXYd,2),RXY(I,3),M(I.l),Md,2),Md,3) 
WRITE(12,123) 

1=1 

1 0  READ(7,*)IPC(I),Y1  (1,1  ),Y2(I,1  ),Y3(I,1 ), 

IF(IPC(I).EQ.999)  GOTO  15 
1=1+1 
GOTO  10 
15  LEN=1-1 

170  FORMAT(lX,9(E15.4,lX)) 
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121  FORMAT(5X,I2,5X,1  '.25X;\  -.2X,1  ,,25X,'I  ',2X,1  -,25X/I ' 

•  2X/I  ',25X,1 ') 

122  FORM AT(12X/r‘,lX,3(F7 .2,1  X)/'T,2X/f’,lX,3(F7.2/lX),,T, 

•  2X,f',lX,3(F7.2,lX);T,2X,f‘,lX,3(F7.2.1X);l) 

123  FORMAT(l2X,t’,25X,'J,,2X/L,,25X,'J',2X,l‘,25X,'J' 

•  2X,*f,25X;l) 

124  FORMAT(lX,E15.4) 


Read  the  VELOCUY.DAT  file  and  prepare  for  isogradient  raytracing. 


DZ=25 
1  =  1 

25  READ(2/,END=30)  TEMP,VV(1) 
1  =  1+1 
GOTO  25 
30  CONTINUE 

GG=(W(M)-VV(I-7))/6. 

DO  35  J  =  1,55 

VV(J)=VV(J-1)+CC 
35  CONTINUE 
LL(1)=12.5 
DEPTH!  1)=0 
DO  40 1=2,55 

LL(I)=LL(I-1)+DZ 
DEPTH(I)=DEPTH(I-1)+DZ 
40  CONTINUE 
DO  43 1=2,55 

G(I-1)=(VV(I)-W(I-1))/DZ 
43  CONTINUE 

C(55)=G(54)+GG/DZ 
DO  45 1=1,52 
45  CONTINUE 
DO80J=U 
N=15 
IFLAG=0 
DO70IT=l,LEN 


Set  variables  the  the  call  to  TGEN  subroutine.  TGEN  returns  the  time  and  elevation. 


A2=Y3(LEN+1,J) 

P1=SQRT((YKIT,J)-Y1(LEN+1,J))**2+(Y2(IT,J)-Y2(LEN+1,J))**2) 

P2=-Y3(IT.]) 

IF(IT.NE.l)  THETA(IT,J)=THETA(IT-1,J) 

CALL  TGEN(LL,G,VV,A2,Pl,P2,THETA(IT,J),TTIME(IT,J), 
V0,V1,DEPTH,IFLAG) 

IFLAG=1 


Find  the  azimuth  angles  from  each  array. 


PH1(JT,J)=DASIN((Y2(IT,J)-Y2(LEN+1.J))/SQRT((Y2(IT,J)- 
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Y2(LEN+1J))**2  +  (YMTJ)-Y1(LEN+1J))**2)) 
IF((Y1  (ITJVYl  (LEN+1  ,J)).LT.O)  PH1(ITJ)=3.141592653S9 
-PHKrrj) 


Locate  the  layer  containing  the  source  and  set  it's  velocity. 


50  IF(-Y30TJ).LE.DEPTH(N).AND.-Y3(rrj).GT.DEPTH(N-l))  GOTO  60 

IF(-Y3(IT,J).GT.DEPTH(N))  THEN 
*  N=N+1 

ELSE 
N=N-1 

JF(N.LE.l)  THEN 
N=1 

GOTO  60 
ENDIF 
ENDIF 
GOTO  50 

60  V(1T,J)=VV(N) 


Calculate  the  BCD  adjust, ,,ents  (Sperical  Coordinates  :  Equtaion  (5)) 


Bl(lT,J)=V(IT,])*DCOS(THETA(rT',J))*DCOS(PHI(IT,J)) 
B2(rr/J)=V(rr,J)»DCOS(THETAaT,J)),DSIN(PHHIT<J)) 
B3(IT  ,])= V  (IT,J)*DSIN(THET  A(ITJ)) 

70  CONT1NUF 
80  CONTINUE 
DO  90  I=1,LEN 
90  CONTINUE 


Calculate  DELTA  and  TNOT. 


D2=0 


T1=0 

DEN=0 

DO  200  FT=1,LEN 

Dl  =(((B1  (IT,1  )-Bl  (IT,2))*(Bl(lT,l)-Bl(IT<2)))+ 
((B2(IT,l)-B2nT,2))*(B2(IT,l)-B2(IT2)))+ 
((B3(IT.l  )-B3(IT,2))*(B3(IT,l  )-B3(IT,2)))) 
DEN=DEN+Dl 
Tl  =Tl  +(IPC(IT)*D1 ) 

D2=D2+(((B1  (IT,1  )-Bl  (IT,2)  W1  (IT,1  )-Yl  (IT,2)»+ 
((B2(IT,1)-B2(IT/2))*(Y2(IT,1)-Y2(IT,2)))+ 
(B3(IT,l)-B3(IT,2>)*(Y3(IT,l)-Y3(rT^)))) 


c(rr,i)  =  BKiT,2)-Bi(rr,i) 
C(IT,2)  =  B2(IT,2)  -  B2(IT,1) 
C(IT,3)  =  B3UT,2)  -  B3(IT,1) 


2a>  CONTINUE 

TD  =  DSQRT(DEN) 


DEN=DEN/LEN 
TNOT=(Tl/ LEND /DEN 
DELTA=-(D2/LEN)/DEN 


Calculate  m  and  g. 


D1=O.ODO 

D2=O.ODO 

DO210IT=l,LEN 

DD=(((B1(IT,1)-Bl  (IT,2))*(B1  (IT,1  )-Bl  (IT, 2)))+ 

•  (B2(IT,1  )-B2(IT,2))*(B2(IT,l  )-B2(IT,2)))+ 

•  ((B3(fT,l  ;-B3(IT,2))*(B3(IT,l  )-B3(IT,2))))* 

*  (IPC(IT)-TNOT)*(IPC(IT)-TNOT)) 

D1  =  D1  +  DD 

D2=D2+(((B1(IT,1  )-Bl  (IT,2))‘(Y1  (!T,1  )-Y10T,2))>-t- 

*  ((B2(IT,1)-B2(IT/2))*(Y2(JT,1)-Y2(IT,2)))+ 

*  (B3(IT,1)-B3(IT,2))*(Y3(1T,1  )-Y3(IT,2))))* 

•  (IPC(IT)-TNOT) 

210  CONTINUE 

MEQ=-(D2/LEN)/(Dl  /LEN) 

TM  =  D1 
DO230J=l,2 

DO  220  IT=1,LEN 

YC1  (rr,J)= Y1  (rrjl-KBl  dT,J)*(DELT  A+(MEQ*(IPC(ID-TNOT)))) 
YC2(IT,J)=Y2(IT,J)+(B2(rrj)*(DELTA+(MEQ*(IPC(IT)-TNOT)))) 
YC3(IT,J)=Y3(IT,J)+(B3(ITJ)*(DELTA+(MEQ*(lPC(rT)-TNOT)))) 
GEQ(IT)=DELTA+(MEQ*(IPC(IT)-TNOT)) 

220  CONTINUE 
230  CONTINUE 


Create  the  output  files. 


SSR  =  0.0D0 
SSB  =  0.0D0 
DO  231  IT  =  1,LEN 

SSR  =  SSR  +  (YC1(IT,1)  -  YC1(IT,2))'*2  *  (YC2(IT,1)  - 

•  YC2UT,2))~2  *  (YC3(IT,1)  -  YC3UT,2))**2 
D1=(((B1(IT,1  )-Bl(IT,2))*(Bl(IT,l)-Bl(IT,2)))+ 

•  ((B2(JT,1)-B2(IT,2))‘(B2(JT,1>-B2(IT,2)))+ 

•  ((B3(IT,1  )-B3(IT,2))*(B3(IT,l  )-B3(IT,2)))) 

SSB  =  SSB  +  D1*CEQ(IT)**2 

231  CONTINUE 

MSR  =  SSR/(LEN  -  2) 

MSB  =  SSB/ 2 
SIC  =  DSQRT(MSR) 

TD  =  TD*DELTA/SIC 
TM  =  MEQ*DSQRT(TM)/S1C 


Calculate  VDEL,SDEL,VM,SDVM  (  variances  and  standard  deviations  ). 
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DUM2  =  O.ODO 
DUM3  =  O.OEX) 

DUM5  =  O.ODO 
DUM6  =  O.ODO 
DO  235 1=1,3 
DO  233  J  =  1,3 
DUM1  =  O.ODO 
DUM4  =  O.ODO 
DO  232  IT  =  1,LEN 

DUM1  =  DUMl  +  (CflTjmiTJ)) 

DUM4  =  DUM4+((IPC(IT)-TNOT)**2  *C(IT,1)*C(IT,))) 

232  CONTINUE 

DUM2  =  DUM2  +  (M(I,J)  *  DUM1) 

DUM5  =  DUM5  +  vM(IJ)  *  DUM4) 

233  CONTINUE 

DO  234  IT  =  1,LEN 

DUM3  =  DUM3  +  (C(IT,1)**2) 

DUM6  =  DUM6  +  ((IPC(IT)  -  TNOT)**2  *  C(IT,I)**2) 

234  CONTINUE 

235  CONTINUE 

VDEL  =  DUM2/DUM3»*2 
SDEL  =  DSQRT(VDEL) 

VM  =  DUM5/DUM6*»2 
SDVM  =  DSQRT(VM) 

DUM2  =  O.ODO 
DO  238  IT  =  1,LEN 
DUM1  =  O.ODO 
EXD  237 1  =  1,3 
DO  236  J  =  1,3 

DUM1  =  DUM1  +  (C(IT,I)*C(IT,J)*M(I,])) 

236  CONTINUE 

237  CONTINUE 

DUM2  =  DUM2  +  ((IPC(IT)  -  TNOT)  •  DUMl) 

238  CONTINUE 

COVDM  =  DUM2/IDUM3  *  DUM6) 

270  READd  1,'(A)',END=280)LINE 
GOTO  270 

280  BACKSPACE  11 

IF(SDEL.EQ.O.ODO)  THEN 
DELTAN  =  SDEL 
MEQN  =  SDEL 
R  =  SDEL 
GOTO  285 
ENDIF 

DELTAN  =  DABS(DELTA)/SDEL 
MEQN  =  DABS(MEQ)/SDVM 
R  =  COVDM /(SDEL*SDVM) 

2 85  DELTA  =  DELTA  *  1000.0DO 
MEQ=MEQ  *1000.000 
SDEL  =  SDEL  *  lOOO.ODO 
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SDVM  =  SDVM  *  1000.0D0 

WRITE(11,350)DELTA,SDEL,DELTAN<MEQ(SDVM<MEQN/R,TNOT(LEN 

3S0  FORMAT(F5.2/2X,F5.4/2X/F6.2(2X/F7.5/2X/F7.6^X/2(F72,2X), 

*  F8.2,2X,I2) 

WRITE(*,*)'  PROGRAM  COMPLETED !' 

END 


SUBROUTINE  TGEN(LL,G,VV,A2,Pl,P2,ANGLE,TIME, 
_  V0,V1,DEPTH,1FLAC) 


DIMENSION  L(55),G(55),V(55),LL(55),W(55) 

DIMENSION  TH(55),T(S5),VZ(55) 

DIMENSION  C2(55),TT(55),DEPTH(55) 

REAL*8  L,G,LL,VV,TH,VZ,C2.TT,T,R0,A1 
REAL*8  A2,P1,P2,C1,C22,THETA,VM 

REAL‘8  THETAZ,FV,R,TIME,ANGLE,EP,DZ,DEPTH,V0,V1,V 


Initialization:  Set  the  value  for  DZ,  the  layer  thickness.  The  sensor  is  assumed  to  be  at 
RANGE  0.  Determine  the  values  for  J,  which  is  1+Number  of  layers  less  than  or  equal  to  the 
sensor  depth,  and  I,  which  is  the  number  of  layers  less  than  or  equal  to  the  torpedo  depth 
Redefine  the  endpoints  of  those  layers  locally  to  be  the  depths  of  the  torpedo  and  sensor. 
Define  local  values  for  the  LL  and  W  arrays. 
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EP=0.1D-6 
DZ=25.0D0 
A1  =0.0D0 
J=1 
1=0 

DO  10  K=l,55 
UK)=LUK) 

V(K)=VV(JO 

IF(DEFrH(K).LE.A2)  J=J+1 
IF(DEPTH(K).LE.P2)  1=1+ 1 
10  CONTINUE 
IFfl.LE.O)  1=1 
N=l+J-I 

V(I)=V(I)+G(I)*(P2-L(I)) 

V(J)=V(J-1)+G(J-1)*(A2-L(J-1)) 

L(I)=P2 

IF(A2.CT.L(J-1))C(J)=G0-1) 

L(J)=A2 

R0=P1-A1 

V1=(V(J)-V(I))/(A2-P2) 
V0=V(J)  -  V1*A2 


Calculate  an  initial  estimate  for  the  angle  and  time  using  a  single  layer  approximation. 


C22=-V0/Vl 
Cl  =«0.5D0)*(P1  -  A1 )) 

C1=C1+((0.5D0)*(L(I)-L(J))*(L(1)+L(J)-2.0D0* 
•  C22)/(P1-A1)) 

IF(IFLAG.GT.O)  C010  48 
THETA=DATAN((A1-C1)/(A2-C22)) 
WRITEC/)'  DEFINE  THETA  AGAIN' 
GOTO  49 

48  THETA=ANGLE 

49  CONTINUE 


Use  the  angle  THETA  to  raytrace  back  through  all  the  layers.  First,  use  the  ray  invariant  (RV) 
and  the  velocity  to  calculate  the  entrance  andgle  at  each  layer.. 


R=A1 


Find  the  maximum  Ss.und  speed. 


VM  =  V(l) 

DO  495  K  =  IJ 

IF  (V(K).CT.VM)  VM  =  V(K) 
495  CONTINUE 

50  RV=DCOS(THETA)/V(J) 
DO60K=IJ 

TH(K)=DACOS<RV*V(K)) 
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T(K)=DT  AN(TH(K)) 

VZ(K)=V(K)-L(K)*G(K) 

C2(K)=-VZ(K)/G(K) 

60  CONTINUE 


Using  the  angle  just  calculated,  iterate  backwords  through  the  layers  from  sensor  to  target  to  get 
the  horizontal  range.  Stop  at  the  depth  of  the  target. 


R=0.0D0 
DO  70  K=J,I+1,-1 

Cl  =R-T  (K)*(UK)-C2(K-1 )) 
R=Cl  +T(K-1  )*(L(K-1  )-C2(K-l)) 
70  CONTINUE 


Test  if  the  value  for  the  range  is  within  torerance.  If  not,  redefine  THETA,  the  initial  angle, 
and  raytrace  again.  If  wthin  tolerance,  calculate  the  time  of  travel  based  on  THETA,  and  return. 


EP=0.1D-6 

IF((DABS(R-Pl)).LE.EP)  GOTO  100 
THETAZ=THETA 

THETA=DATAN(DTAN(THETAZ)*(R-A1)/R0) 

GOTO  50 

100  TT(J)=DLOG((l  .0D0+DSIN(TH(J)))/(DCOS(TH(J)))) 
TIME=0.0D0 
DO  110  K=J-1,I,-1 

TT(K)=DLOG((1.0D0+DSIN(TH(K)))/(DCOS(TH(K)))) 

TIME=TIME+CTT(K)-TT(K+1))/G(K) 

110  CONTINUE 

ANCLE=THETA 

RETURN 

END 
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